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EXECUTIVE  SUMMARY 


This  report  describes  theoretical  woik  carried  out  in  support  of  the  Plasma  Radiation  Source 

(PRS)  DECADE  Program  by  Code  6720  in  FY  1993.  The  six  sections  of  the  rqwit  discuss  the 

following  subjects: 

(1)  In  section  I,  DECADE’S  predicted  PRS  performance  is  discussed.  During  FY  93,  a  capability 
to  evaluate  the  effect  of  different  DECADE  front-end  designs  on  Marx-bank-energy-to-load- 
Idnetic-energy  conversion  efficiency  was  developed,  and  a  significant  advancement  was  made 
in  our  ability  to  predict  how  load  kinetic  energy  conversion  to  K-shell  x  rays  scales  with  mass 
and  implosion  velocity.  This  advancement  was  made  possible  by  a  systematic  comparison  of 
theoretically  predicted  K-shell  yields  with  measured  yields  in  experiments  on  Double  Eagle, 
Saturn,  and  Phoenix  machines.  The  lattn  expoiments  were  designed  this  year  e:q)ressly  for 
this  task.  There  are  two  ways  to  qrproach  the  DECADE  coupling  and  scaling  problems.  In 
one  iqjproach,  emphasis  is  placed  on  the  technical  issues  associated  with  power  flow  to  the 
PRS  load  and  on  the  design  of  multi-megajoule  pulse-power  generators  for  use  as  PRS  x-ray 
simulators.  In  the  other  ^roach,  emphasis  is  placed  on  the  load-specific  requirements  for 
efficiently  converting  electrical  energy  into  x  rays  (e^cially  x  rays  that  are  tailored  for  use  in 
NWET).  The  first  qjproach  to  DECADE  guides  the  design  of  DECADE’S  PRS  front-end.  In 
this  case,  it  is  assumed  that  if  a  pulse  power  machine  can  be  built  to  meet  the  power  flow  and 
energy  coupling  PRS  requirements,  then  x  rays  will  be  produced  as  predicted  by  the  scaling 
laws.  This  assumption  allows  emphasis  to  be  placed  on  machine  design  problems.  In  FY93, 
these  scaling  laws  were  put  on  a  more  solid  foundation.  The  culmination  of  several  recent 
years  of  experimental  and  theoretical  efforts  to  improve  the  scaling  laws  is  summarized  in 
Section  I.  As  a  consequence,  x-ray  yield  scaling  in  aluminum  can  now  be  predicted  over  three 
or  more  orders  of  magnitude,  and  machine  requirements  can  be  defined  for  achieving  greater 
simulation  fidelity.  On  scales  of  several  orders  of  magnitude,  a  factor  of  two  difference  in 
yields  can  be  ignored.  However,  when  a  DECADE-  or  Jupiter-dass  PRS  machine  is  built,  a 
factor  of  two  difference  in  performance  will  be  very  significant  in  terms  of  doUars^’hoton. 
Thus,  section  I  concludes  with  a  discussion  of  some  non-machine  Goad  oriented)  problems 
that  were  addressed  as  part  of  the  FY93  technical  risk  reduction  effort  for  DECADE.  In  both 
cost^oton  and  simulation  fidelity  areas,  the  load  physics  PRS  issues  become  important. 

(2)  The  scaling  relations  being  used  to  predict  DECADE’S  PRS  performance  capabilities  need 
theoretical  and  experimental  validation  for  moderate  atomic  number  loads,  (18  <  Z  <  29). 
One  way  to  validate  theoretically  is  to  carry  out  a  set  of  calculations  for  titanium  or 
copper  that  is  similar  to  those  that  were  done  for  aluminum.  However,  experimental 
implosions  are  much  softer  G*e.,  much  lower  in  density  and  much  fatter  in  extent)  than  the 
aluminum  implosions  that  were  calculated  in  deriving  the  earliest  version  of  the  K-shell 
scaling  laws.  These  calculations  were  carried  out  without  the  use  of  phenomenologically 
determined  enhancements  to  classically  calculated  plasma  viscosities,  heat  conductivities, 
and  electrical  resistivities.  During  FY93,  a  significant  effort  was  completed  to  determine 
dwse  enhancements  from  PI  experimental  data  in  order  to  obtain  a  more  realistic  scaling  of 
K-shell  yield  with  load  mass  and  implosion  velocity.  In  particular,  the  shift  in  brealqwint 
mass  was  determined  from  an  analysis  of  argon  and  aluminum  experimental  data.  This 
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wodc  is  discussed  in  Section  U.  Because  soft  implosion  modeling  generally  produces  better 
agreement  between  theory  and  experiment,  it  also  creates  new  problems  of  discriminating 
between  physical  situations  where  conventional  slug  modeling  and  current-off  procedures  are 
valid  and  where  they  are  not.  These  problems,  which  are  the  subject  of  current  investigation, 
are  also  discussed  in  Section  n.  They  are  relevant  to  the  problem,  also  under  investigation,  of 
increasing  PRS  emissions  in  the  10  to  20  ImV  region  on  DECADE  or  ACE  IV. 

(3)  IWo  problems  arise  as  machine  energy  is  scaled  upwards.  On  the  one  hand,  the  PRS  loads 
must  be  imploded  from  larger  distances,  leading  to  a  possible  degradation  in  implosion 
symmetry  due  to  the  onset  of  Rayleigh-Taylor  instabilities.  On  the  other  hand,  larger  amounts 
of  mass  must  be  imploded,  leading  to  questions  of  how  nozzle  designs  should  be  scaled.  The 
Implication  of  a  2-D  MHD  fluid  code  to  address  these  problems  is  discussed  in  Section  m.  It 
is  found,  for  example,  that  the  growth  of  plasma  instabilities  can  be  mitigated  by  distributing 
a  gas  PRS  load  over  over  a  large  volume  and  imploding  it  more  as  a  snowplow  than  a  plasma 
shell.  However,  this  option  leads  to  less  kimtic  energy  generation  than  a  shell  implosion. 
The  tradeoff  between  instability  reduction  and  maximizing  kinetic  energy  is  an  important  2-D 
issue  under  investigation.  The  2-D  code  is  also  being  used  to  investigate  the  influence  of 
different  gas  puff  nozzle  designs  on  PRS  performance.  This  gas  puff  work  is  being  done  in 
support  of  PI  Double  Eagle  experiments,  which  provides  a  benchmark  before  they  are  ipplied 
to  DECADE. 

(4)  The  POS  is  an  important  circuit  element  in  DH14DE  when  it  is  operated  in  bremsstrahlung 
mode.  Its  role  in  PRS  mode  is  less  certain.  Over  the  last  coipleofyears,  an  NRL  transmission 
line  model  of  an  inductive  energy  store  (lES)  generator  was  constructed.  Its  POS  model  was 
benchmaiked  in  FY92  to  NRL  Hawk  POS  data  and  qplied  in  FY93  to  the  study  of  power 
flow  in  DECADE.  In  Section  IV,  the  DECADE  transmission  line  model  is  presented,  and 
results  of  its  benchmarking  against  Pi’s  D^ADE  model  are  discussed.  When  coupled  to 
a  gas  bag  model  PRS  load,  the  NRL  model  supports  flie  idea  that  the  DECADE  generator 
should  be  able  to  meet  DECADE’S  radiation  specs  in  PRS  mode  with  or  without  a  POS. 
With  a  POS,  however,  higher  yields  should  be  achievable  with  copper  or  krypton  loads  than 
without  provided  the  POS  operates  at  a  sustainable  high  quality  level.  This  prediction  was 
tested  theoretically,  in  part,  by  a  sensitivity  study  of  energy  transfer  to  a  PRS  load  to  variations 
in  the  POS  performance  parameters. 

(5)  Predictions  of  DECADE’S  performance  and  of  load  designs  to  improve  its  performance  rely 
on  the  ability  of  theoretical  models  to  accurately  describe  PRS  load  dynamics  and  on  the 
ability  of  PRS  experiments  to  to  meet  performance  objectives  by  design.  To  increase  the 
reliability  of  theory  and  experiment,  procedures  to  bmchmatk  z-pinch  calculations  against 
erqrerimental  data  must  be  developed.  Some  efforts  in  this  direction  are  described  in  Section 
V.  A  set  of  e]q)aiments  was  recently  carried  out  at  flie  Weizmann  Institute  of  Science  in 
Israel  in  which  an  extensive  set  of  measurements  were  made  of  the  time-dependent  radial 
velocity  distributions  of  sin^y  to  five  times  ionized  ions  in  an  imploding  plasma  shell.  It 
was  speculated  that  an  ionization  wave  propagating  much  faster  than  the  local  radial  ion 
velocities  was  observed.  In  Section  V,  an  analyas  is  described  in  which  data  similar  to  that 
observed  in  the  Israli  e3q)miments  is  obtained  by  post  analyzing  a  1-D  MHD  calculation. 


The  correspondence  of  this  data  to  the  Israli  data  is  discussed.  The  development  of  similar 
diagnostic  procedures  is  planned  for  FY94  to  be  used  in  the  analysis  of  the  Saturn,  PI,  and 
Maxwell  ACE  IV  data. 

(6)  Finally,  a  small  effort  continued  in  FY93  to  construct  scalable  L-shell  models  for  use  in 
developing  diagnostics  for  PRS  experiments  on  DECADE  and  ACE  IV.  This  work  is  also 
relevant  to  problems  such  as  L-shell  bum-through.  The  model  development,  which  is 
discussed  in  Section  VI,  is  patterned  after  the  successful  use  of  hydrogen-  to  helium-like 
line  ratios  to  derive  information  about  the  state  of  a  plasma  radiating  in  the  K-shell.  The 
generalization  of  these  K-shell  diagnostic  procedures  to  the  L-shell  is  accomplished  by  a 
substitution  of  the  closed  shell  neon-like  ionization  stage  for  the  closed  shell  helium-like 
stage  and  the  fluorine-like  for  hydrogen-like  stage.  In  the  L-shell,  the  multiplet  structure 
of  the  bound  states  complicates  the  problem  of  building  reliable  L-shell  dynamical  models; 
however,  resources  are  available  with  which  to  address  these  problems.  The  fluorine-like 
model  development  described  in  Section  VI  makes  use  of  atomic  code  capabilities  were 
developed  at  die  Los  Alamos  National  Laboratory. 
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I.  PRS  ON  DECADE 


A.  DECADE  ENERGY  CONSIDERAnONS 

Designing  DECADE  for  use  as  a  plasma  radiation  source  (PRS)  involves  two  problems:  (1) 
the  efficiency  of  coiq)ling  Marx  bank  enogy  to  the  PRS  load  must  be  realistically  evaluated  for 
different  DECADE  PRS  front-end  designs,  and  (2)  the  efficiency  of  converting  coupled  electrical 
energy  into  K-  and  L-shell  kilovolt  radiation  must  be  realistically  determined  for  a  variety  of  load 
designs  and  for  a  selection  of  low,  moderate,  and  hi^  Z  elements.  During  FY  93,  a  usability  to 
evaluate  die  effect  of  different  front-end  designs  on  DECADE  Marx  bank/PRS  coupling  efficiency 
was  devdoped,  and  a  significant  advancement  was  made  m  our  understanding  of  the  scaling 
of  K-shell  yields  with  mass  and  implosion  velocity  to  DECADE.  These  developments  will  be 
described  in  diis  section.  The  latter  advancement  was  made  possible  by  a  systematic  comparison  of 
theoretically  predicted  K-shell  yields  with  measured  yields  on  Double  Eagle,  Saturn,  and  Phoenix 
machines. 

When  emphasis  is  placed  on  technical  issues  associated  with  the  design  of  multi-megajoule 
pulse-power  generators,  a  machine  rqiproach  to  DECADE  PRS  performance  is  needed  to  guide  the 
design  of  the  DECADE  generator’s  PRS  front-end  (with  or  without  a  POS).  A  PRS  load  can  then 
be  defined  in  simplest  terms  by  the  amount  of  mass,  m,  to  be  imploded  and  by  the  initial  radius,  , 
from  which  it  is  imploded.  Both  of  these  quantities  are  assumed  to  be  completely  free  variables. 
The  most  basic  machine  design  problem  concerns  the  determination  of  the  amount  of  useful  energy 
that  the  machine  can  deliver  to  the  PRS  load  for  omversion  into  x  rays.  The  simplest  way  to 
address  tiiis  problem  is  to  use  a  (lumped)  circuit  model  that  describes  a  given  DECADE  machine 
to  couple  to  the  PRS  load.  A  slug  implosion  model^  of  the  load  can  then  be  used  to  calculate  the 
kinetic  energy  imparted  to  the  given  load  by  the  giv«a  machine.  These  calculations  are  stopped  at 
a  pre-determined  radius  (2  nun)  and  the  firuil  imploaon  velocity,  Vtmp,  of  the  load  is  determined. 
The  choice  of  2  mm  is  governed  by  two  factors:  (1)  the  further  in  the  slug  implosion  proceeds, 
tite  larger  Vimp  becomes  and  (2)  the  further  out  the  implosion  is  stopped,  the  more  likely  it  is  tiiat 
plasma  back-pressure  carmot  halt  the  acceleration.  Finally,  the  results  of  these  calculations  are 
displayed  on  an  m-v^mp  contour  plot. 

Fig.  (1)  shows  such  a  plot.  The  darkened  crescent  region  in  Fig.  (1)  represents  the  ^an  of 
implosion  velocities  that  were  calculated  using  the  circuit  model  for  DECADE  that  is  shown  in 
Fig.  (2)  for  the  range  of  radii  and  masses  shown:  2  cm  <  r,-  <  3.S  cm  and  0.22  mg/cm  <  m  < 
10.6  mg/cm.  Diis  darkened  region  overlqjs  witii  three  regions  that  are  labeled  in  Fig.  (1)  as 
regions  of  inefficient,  efficient,  and  inefficient  x-ray  production  respectively.  The  labeled  regions 
in  this  figure  are  for  titanium.  Because  the  darkmed  crescent  overl{q>s  the  efficient  emission 
region,  titanium  loads  ffiould  be  able  to  efficiently  convert  implosion  kinetic  energy  to  K-shell  x 
rays  on  DECADE,  but  oidy  for  a  certain  choice  of  radii  and  masses.  More  information  about  the 
origin  of  tiiis  figure  will  be  published  in  an  NRL  memorandum  report  entitied  “PRS  Projections 
for  DECADE”. 

The  model  in  Fig.  (2)  contains  the  following  basic  elements  of  a  DECADE  circuit:  the  Marx 
goierator,  Ci,  the  transfer  ciqiacitor,  Cj,  the  transfer  switch,  iZs,  the  pulse  forming  line,  Ls,  the 
plasma  opening  switch  (PEGS),  Rpeot,  the  magnetically  insulated  transmission  line  (MTTL)  of  the 
nUS  diode,  £4 ,  and  tiie  PRS  load,  Viood  •  These  elements  were  given  the  values  shown  in  the  figure. 
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From  Aus  lumped  circuit  of  DECADE,  one  can  compute  the  efficiency  of  converting  the  initially 
stored  energy  of  the  Marx  bank,  EMan,  into  PRS  kinetic  energy.  For  example,  for  the  same 
range  of  masses  and  load  radii  that  were  used  to  obtain  the  span  of  implosion  velocities  shown  in 
Fig.  (1),  one  obtains  the  corresponding  span  of  Maix-bank-to-load  conversion  efficiencies  shown 
in  Fig.  (3)  (for  a  3  cm  long  array).  For  the  Fig.  (2)  circuit,  the  largest  mass  loads,  m  ~  10  mg/cm, 
are  also  the  most  efficient  converters  of  Marx  energy  into  load  kinetic  energy  ('^  1S%*17%).  But, 
as  Fig.  (2)  shows,  this  mass  of  titanium  on  DECADE  is  an  inefficirat  K-shell  radiator.  Efficient 
emission  is  predicted  to  occur  only  for  0.6  mg/cm  <  m  <  1.8  mg/cm.  For  these  masses,  the 
conversion  efficiency  of  Marx  bank  energy  to  kinetic  orergy  is  betweoi  8%  and  15%. 

One  can  anticipate  that,  at  the  present  time  under  optimum  conditions,  30  to  50%  of  the  kinetic 
energy  generated  during  a  PRS  implosion  can  be  converted  into  K-shell  x  rays  in  the  efficient 
x-ray  production  region  of  m  -  Vimp  space.  By  multiplying  together  the  above  two  conversion 
efficiencies,  one  concludes  that  2.4  to  7.5%  of  the  Marx  bank  energy  can  be  converted  to  K-shell 
X  rays  by  the  Fig.  (2)  circuit  model  (8%  to  15%  x  30%  to  50%).  Since  Euarz  =  8.5  MJ,  these 
considerations  predict  that  the  Fig.  (2)  DECADE  machine  should  be  able  to  generate  between  200 
to  640  kJ  of  titanium  K-shell  radiation.  This  factor  of  three  in  x-ray  output  could  be  extremely 
I  important  to  the  simulator  users  community. 

B.  K-SHELL  SCALING,  THEORY  VERSUS  EXPERIMENT 

IWo  implosion  kinetic  energies  determine  the  radiation  characteristics  of  PRS  loads: 

(1)  the  kinetic  energy  per  ion,  Ki  =  (l/2)mivf^p,  and  the  kinetic  energy  per  centimeter, 

Ke  =  where  mi  is  the  mass  of  an  ion.  Ki  determines  the  ability  of  the  load  to  radiate 

from  the  K-shell.  L-shell  radiation  losses  are  ignored,  then  Ki  must  be  larger  than  a  predetermined 
minimum  energy,  Emin,  which  one  can  estimate  by  the  formula,  Emin  =  1.012Z®  *®*,  where  Z 
is  the  atomic  number  of  the  PRS  load.  Since  K-:diell  radiation  is  also  emitted  per  ion  in  units  of 
Emin  in  the  efficient  emission  region,  it  is  convenient  to  define  a  dimensionless  kinetic  energy  per 
ion,  Tj,  defined  by  17  =  Ki/Emin-  One  can  then  work  interchangeably  in  terms  of  the  variables 
{m,vimp}  or{m,i7}. 

Our  present  understanding  of  K-shell  yield  scaling  evolved  over  a  period  of  the  last  five  years 
from  a  concerted  theoretical  and  experimental  effort  between  NRL  and  Physics  International,  Inc.. 
Roughly  speaking,  it  involved  the  following  four  separate  steps:  (1)  A  simple  slug  model  theory 
of  K-shell  yield  scaling  was  developed,  and  used  to  plan  a  set  of  1-D  radiative  hydrodynamics 
implosion  calculations  that  systematically  quantify  the  radiative  capabilities  of  z-pinches.  In  these 
calculations,  m  was  held  constant  and  -q  was  varied  or  17  was  held  constant  and  m  was  varied. 

(2)  The  results  from  these  calculations,  when  plotted  on  a  log-log  graph,  showed  there  were 
roughly  two  scaling  regimes  of  K-shell  yield  with  mass,  which  made  sharp  transitions  into  one 
another.  The  machine  requirements  for  producing  K-shell  radiation  from  elements  having  atomic 
numbers  higher  than  aluminum  were  then  predicted.  (3)  Detailed  comparisons  were  next  made 
of  the  calculated  z-pinch  results  with  carefully  analyzed  experimental  spectral  data.^  (4)  The 
calculations  were  then  systematically  modified  this  year  by  enhancing  the  plasma  viscosity,  the 
heat  conductivity,  and  the  electrical  resistivity  through  the  use  of  multiplicative  factors  in  order  to 
improve  their  agreement  with  PI  e:q)eriments.  They  had  the  effect  of  significantly  modifying  the 
scaling  predictions.’  The  newly  modified  scaling  relation  was  used  this  year  to  analyze  two  sets 
of  recently  completed  Saturn  and  Phoenix  aluminum  experiments.  The  Phoenix  experiments  were 
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designed  and  analyzed  this  year.  Both  the  Saturn  and  Phoenix  experiments  were  ttesigned  to  probe 
K-shell  emission  through  the  predicted  efficient  emission  region  and  to  compare  to  the  theory  that 
was  developed  firom  tlw  analysis  of  Double  Eagle  experimmits.  The  analysis  of  these  experiments 
led  to  the  (tevelopment  of  a  new  scaling  relation  with  which  to  predict  DECADE  performance. 

Tlw  expmiments  at  Sandia  National  Laboratories  and  at  NSWC  on  the  Saturn  and  I%oenix 
machines  respectivdy,  make  it  possible  to  test  the  yield  predictions  over  two  different  ^ans  of  m 
and  Vimp  in  m  -  vimp  space.  Hie  Saturn  experiments  were  performed  at  the  end  of  FY92.  The 
hi^iest  yielding  of  the  different  mass  load  Saturn  experiments  are  listed  below  in  a  table.  Each 
array  had  a  length  of  2  cm: 


SATURN  PRS  Experiments 


Wire  Diam 
(ndl) 

Wire  Mass 
(pg/cm) 

Array  Diam 
(cm) 

H 

Wires 

# 

K.E. 

kJ 

Yield 

kJ 

C.E. 

% 

0.6 

120 

3.04 

22 

24 

230 

8.0 

3.5 

0.7 

160 

2.81 

17 

24 

230 

12.2 

5.3 

0.75 

185 

2.73 

15 

24 

230 

22.4 

EBI 

1.0 

330 

1.72 

6.0 

24 

170 

76.0 

44.7 

1.2 

470 

2.29 

5.6 

24 

230 

61.2 

26.6 

1.5 

740 

1.46 

EQI 

24 

170 

30.0 

17.6 

1.7 

950 

1.4 

22 

24 

170 

20.0 

11.8 

imiQQiiniii 

1310 

1.34 

1.6 

24 

170 

14.2 

8.4 

2J 

2050 

125 

0.95 

24 

170 

4.9 

2.9 

In  each  of  these  experiments,  24  two  centimeter  long  wires  of  different,  but  standard,  wire 
diameters  were  mounted  at  the  different  array  diameters  listed  in  the  table.  The  total  kinetic 
energies  and  the  17  values  diat  are  listed  were  computed  from  a  slug  model  using  the  Saturn 
lumped  circuit  as  its  driver.  The  (K-shell)  total  yields  were  measured,  and  each  percent  conversion 
efficiency  (C.E.)  diat  is  listed  in  die  table  is  the  rado  of  measured  yield  to  computed  kinedc  energy 
in  each  experiment.  For  the  similar  set  of  experiments  carried  out  at  Phoenix  in  FY93,  again,  only 
the  highest  yielding  shots  are  listed. 

The  specific  locadon  in  m  -  Vimp  space  of  these  experiments  is  shown  in  Fig.  (4)  along  with 
contours  for  predicted  K-shell  yields  Hie  nine  dots  represent  the  nine  2  cm  load  Saturn  experiments 
and  dw  nine  triangles,  the  nine  2  cm  load  Phoenix  experiments.  The  K-shell  yields  that  were 
measured  in  each  of  these  experiments  are  plotted  in  Figs.  (5)  and  (6)  respeedvely.  In  both  sets 
of  experiments,  the  maximum  yield  was  recorded  in  the  eiqieriment  that  lay  within  the  efficient 
region  but  was  closest  to  the  m  =  rngp  boundary  curve.  Moreover,  as  predicted,  the  yields  fell  as 
the  experiments  generated  more  implosion  velocity,  but  moved  further  from  the  efficient  region. 
Similarly,  as  the  implosicm  velocity  was  decreased  and  the  experiments  moved  in  the  direcdon  of 
die  17  =  1  lower  boundary  line,  the  yields  fell.  The  yield  contours  that  are  drawn  in  Fig.  (4)  were 
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array  length  s  2  cm 


Phoenix  Aluminum  Experiments 


Wire  Diam 
(mil) 

Wire  Mass 
(/ig/cm) 

Array  Diam 
(cm) 

H 

Wires 

# 

K.E. 

kJ 

Yield 

kJ 

C.E. 

% 

0.6 

79 

2.62 

7.0 

16 

46 

8.6 

18.7 

0.7 

107 

23% 

Bl 

16 

52 

19.5 

313 

0.8 

140 

225 

4.05 

16 

48.5 

26 

54 

1.2 

158 

2.35 

3.8 

8 

52 

21.7 

41.7 

1.0 

219 

2.1 

2.6 

16 

55 

IS 

33 

13 

246 

1.9 

2.2 

8 

46 

is 

49.8 

1.6 

280 

2.0 

2.0 

8 

49 

21 

43 

1.2 

315 

2.26 

2.2 

16 

60 

14.7 

24.5 

1.6 

560 

1.15 

0.6 

16 

31 

6.5 

21 

derived  £nmi  hard  implosion  calcularions  (corrected  only  for  the  mass  brealqxnnt  shift)  and  do  not 
show  dw  drop  in  yield  with  rj  that  is  seen  in  the  data.  However,  soft  implosion  calculadms  do 
show  this  trend. 

An  importam  difference  showed  up  between  the  Phoenix  and  the  Saturn  data.  In  all  of  the 
six  Phoenix  shots  that  lay  within  the  efficient  region,  the  percent  conversitm  of  kinetic  energy 
into  K’SheU  x  rays  was  greater  than  die  theoretically  conservative  value  of  30%.  For  Saturn, 
however,  only  the  maximum  yielding  of  the  five  experiments  within  the  efficient  region  exceeded 
30%.  However,  in  this  (me  case,  the  Saturn  conversion  efficimcy  of  45%  was  close  enough  to 
the  Phomix  maximum  ccmversion  efficiency  of  34%  to  lie  within  the  eiqierifflental  variation  of  the 
two  machines. 

A  tentative  conclusicm  can  be  drawn  from  the  location  of  the  maximum  yielding  of  the  Saturn 
and  Phoenix  mcpeiiments.  Since  they  both  lie  close  to  the  curve,  vk  =  vbp  =  yK{‘m)\m=mBP  = 
yBPi^Bp)f  one  can  use  this  curve  to  predict  an  optimal  aluminum  performance  for  DECADE 
given  by 

VBP  =  0.0056m^j,^S  (1) 

when  /k  =  0.5.  This  curve  is  drawn  in  Fig.  (7),  and  the  Phoenix  and  Saturn  data  points  that 
represem  die  maximum  yielding  experiments  are  included  (m  the  figure.  The  DECADE  and 
Jupiter  points  were  placed  cm  the  curve  by  choosing  the  largest  mass  that  can  be  accelerated  to  the 
brealqmint  velcxaty  by  these  machines.  For  DECADE,  this  mass  was  taken  from  Fig.  (1). 

C.  ACE  IV  EXPERIMENTS 

It  is  expected  that  ACE  IV  of  Maxwell  Laboratories  will  come  on  line  within  the  coming 
year.  This  machine  is  a  direct  drive  monolithic  Inductive  Energy  Store  (lES)  genoator.  Since  it 
is  die  first  multi-megamp  current  lES  driver,  it  is  clearly  useful  to  do  PRS  eiqieiiments  <m  ACE 
rV  in  smqxm  of  die  DECADE  PRS  front-end  d^gn.  In  anticipation  of  such  eiqierimaits,  we 


4 


calculated  the  K-shell  yields  from  aigon  puff  gas  loads  on  ACE  IV  over  a  range  of  initial  radii  and 
mass  loadings  using  a  transmissitm  line  model  to  describe  the  loadMachine  coupling.  We  found 
dutt  for  an  ^ipropiiate  mass  loading  and  for  a  pinch  length  of  3  cm,  the  total  argcm  K-shell  yield 
diould  range  from  40  kJ  for  a  initial  radius  of  0.75  cm  to  100  kJ  at  an  initial  radius  of  6  cm.  This 
significant  impcovemoit  in  yield  as  one  moves  out  in  initial  radius  could  be  used  as  a  test  of  large 
radii  implosimis  to  determine  if  any  degradation  occurs  due  to  Rayleigh-Taylor  type  instabilities 
during  run  in. 

The  calculati<m  of  the  K-shell  radiation  yield  is  adqrted  frmn  Whitney,  et  al.  scaling  (J.  .^ipl. 
niys.,  67,  p.l725,  1990).  Recently,  diis  theory  has  beat  reworked  so  that  the  yield  Yk  can  be 
treated  as  a  function  tmly  of  die  final  implosion  velocity  v/,  die  mass  loading  per  unit  length  Mjt, 
and  the  atomic  number  2n«e.  Omtours  of  constant  yield  over  the  v/  -  M/l  plam  are  presented 
for  aigtm  in  Fig.  (8).  For  purposes  of  evaluating  the  predictive  disability  of  such  a  grqdi,  the 
measured  peak  argcm  K-shell  yields  (per  cm)  <m  DOUBLE  EAGLE  and  SATURN  are  noted  as 
circles.  In  die  inefficient  regime,  die  yield  scales  as  (Af/f)^,  while  in  the  efficiem  regime  it  scales 
as  {M/l).  The  K-diell  scaling  theory  becomes  invalid  when  die  kinetic  energy  per  particle  divided 
the  amgy  needed  to  attain  K-shell  ionizatimi  conditions  is  less  than  or  close  to  unity.  This 
.  ratio  is  termed  17  in  the  reference  above,  and  two  values  are  contained  in  the  plot. 

As  presented,  the  K-shell  yield  scaling  is  independent  of  the  driver,  demanding  only  a 
knowledge  of  the  final  velocity  and  mass  loading.  In  order  to  predict  yields  for  ACE  IV,  we  used 
the  transmissitm  line  circuit  model  tor  ACE  IV  pictured  in  Fig.  (9).  It  has  been  compiled  frtnn 
information  forwarded  to  Code  6720  by  E.  Waisman  and  R.  Ingermanson  of  Maxwell  lalxiratories. 
For  the  plasma  opening  switch  (POS)  we  used  their  results,  which  show  a  quadratic  rise  in  POS 
impedance  to  1.3  ohms  in  80  nsec.  If  we  assume  that  the  argon  gas  puff  is  initially  a  thin  shell, 
one  can  drive  a  slug  model  implosimi  with  this  circuit  and  stop  tlw  implosion  once  the  radius  of 
the  slug  reaches  1/10  of  its  initial  radius.  The  dug  velocity  at  this  point  is  equated  with  the  final 
velocity  v/.  For  ACE  TV  we  todc  die  pinch  length  as  3  cm.  The  current  profiles  for  various  values 
of  mass  loading  times  the  square  of  tiie  initial  radius,  Af  A’,  are  contained  in  Fig.  (10).  For  a  dug 
model,  the  current  and  implodcm  time  ate  only  dependoit  upon  this  produa.  Note  tiiat  for  A/i2^  > 
100  mg  cm^ ,  the  peak  current  into  the  load  is  8  MA,  and  that  the  inductive  notch  decreases  slightly 
for  longer  implodon  times. 

Vfith  the  cdculated  fiiud  velocity  and  the  mass  loading,  which  is  input  at  the  beginning  of 
the  calculation,  one  can  employ  Fig.  (8)  to  determine  the  radiation  yields  for  argon.  Contours  of 
constant  K-shell  yield  vuot  yield  per  cm)  ate  presented  in  Fig.  (11)  as  dotted  lines.  This  is  the 
plane  of  initid  radius  Ro  versus  mass  loading  Af ,  and  rq>resents  the  plane  of  initid  conditions 
assuming  a  dug  implodon  modd  is  an  adequate  description.  Lines  of  fixed  q  vdues  are  shown  as 
solid.  The  cmistant  dope  lines  ate  labded  with  their  corre^nding  vdue  of  MR^  in  mg  cm^,  and 
can  be  used  in  ctmjunction  with  Fig.  (10)  to  estimate  implodon  times.  It  is  imperative  to  state  that 
the  stupe  of  the  curves  in  tins  figure  ate  very  depoident  ipon  the  machine  driving  tiie  implodon. 
No  generd  condudons  can  be  made  regarding  otiier  generators. 

To  summarize  this  andysb,  tiie  peak  yields  for  argon  on  ACE  IV  as  configured  in  Fig.  (9) 
tqppear  to  arise  for  17  between  2  and  3.  Condder  a  totd  mass  loading  of  1.5  mgm.  If  tiie  initid 
radius  is  0.75  cm,  one  can  oqiect  ~40  kJ  of  K-shell  radiation  Yk.  On  the  other  hand,  moving 
out  to  6  cm  initid  radius  more  than  doubles  the  Yk  and  increases  the  implodon  time  from  ^80 
to  ~340  nsec.  It  would  be  useful  to  run  a  series  of  mtperiments  of  ACE  TV  where  the  mass  in 
tiw  argon  gas  puff  is  constant  while  the  outer  radius  is  substantially  increased.  Information  on 
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the  stability  of  large  radius,  long  time  implosions  would  have  direct  implications  for  the  design 
cmidititms  and  yield  expectations  of  the  next  generation  of  PRS  drivers. 

D.  PRS  ISSUES 

Load  oriented  issues  must  also  be  addressed  in  order  to  optimize  DECADE’S  performance  as 
an  NWE  simulator.  We  will  briefly  discuss  some  of  these  issues,  and  mention  some  of  the  woric 
dut  was  dcme  this  year  in  addressing  them.  They  are  subsumed  under  the  following  four  topics: 

•  PRS  Initial  Conditicms 

•  PRS  On-axis  Symmetry  and  Implosion  Stability 

•  20  to  30  KeV  Photons 

•  Radiative  ColltQ>se 

Each  of  these  topics  is  discussed  briefly  in  turn. 

PRS  Initial  Conditions 

IWo  basic  assumptions  underlie  the  scaling  calculations  described  above.  All  of  the 
,  hydrodyrutmics  calculations  begin  with  a  given  mass  of  pre-ionized  cylindrically  symmetric 
plasma  ti^tly  located,  on  average,  at  the  initial  position  of  a  wire  array.  They  assume  (1)  that  the 
power  flow  from  the  pulsed  power  machine  can  be  coupled  to  the  wire  load  to  uniformly  ionize  it 
quickly  and  (2)  that  the  implosion  then  takes  place  with  all  of  the  wire  mass  imploding  inwards  as 
described  by  a  dug  model  driven  by  the  j  x  B  forces  of  the  current  discharge.  These  assumptions 
are  not  always  justified.  Large  diameter  wires  do  not  necessarily  ionize  uniformly,  some  of  the  wire 
mass  may  be  left  behind  in  the  implosion,  and  mass  can  be  driven  in  under  a  combination  of  slug 
and  snowplow  dynamics.  The  first  two  of  these  wire  phenomena  are  suggested  by  experiments, 
and  the  fliird,  by  flie  analysis  of  PI  data  that  was  carried  out  this  year.’  One  conclusion  that  can  be 
drawn  from  these  considerations  is  that,  as  machine  power  is  increased,  the  problem  of  evaluating 
the  initial  response  of  a  wire  load  to  the  current  pulse  could  be  the  most  significant  part  of  the 
power  flow  problem. 

Other  evidence,  shown  in  Figs.  (12)  and  (13),  was  found  this  year  to  suggest  the  importance 
of  the  initial  breakdown  response  of  a  PRS  load  to  the  overall  performance  of  an  x-ray  simulator. 
In  Fig.  (12),  the  Saturn  yield  data,  discussed  above,  is  plotted  as  a  function  of  array  diameter.  This 
figure  also  includes  Double  Eagle  data.  The  dip  in  the  Saturn  yield  at  an  array  diameter  of  2.2  cm 
is  correlated  to  the  rise  in  the  wire  i<m  density  generated  by  the  wire  explosion  halfway  into  the 
current  ri%  (Fig.  (13)).  A  self-similar  hydro^mamics  model  was  used  to  calculate  the  erq>losion 
dynamics.  It  assumed  that  the  wire  expansion  was  uniform  and  that  this  expansion  was  driven  by 
a  time  evolving  isothermal  plasma.  This  work  is  described  in  greater  detail  in  Ref.  (4). 

PRS  On-axis  Symmetry  and  Implosion  Stability 

As  machine  power  is  increased,  the  initial  array  diameter  of  a  PRS  load  will  also  need  to  be 
increased  for  sevnal  reasons.  Large  radii  implosions  will  reduce  some  of  the  difficulties  in  getting 
power  to  flow  initially  into  the  diode  and  to  couple  effectively  to  the  PRS.  They  will  rdso  reduce 
the  requirement  for  a  fast  rising  current  pulse.  Finally,  they  should  allow  higher  aspect  ratio  (i.e. 
larger  ratios  of  initial  to  final  array  diameter)  implosions  to  be  achieved.  Increases  in  this  rat  o  will 
trandate  mto  increases  in  Vi^p  and  to  a  greater  likelihood  fliat  the  K-^dl  of  high  Z  elements  can 
be  ignited  by  these  implosions. 
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As  array  diameters  are  increased,  it  could  become  increasingly  difficult  to  maintain  the 
symmetry  of  the  implosion,  which  is  critical  for  achieving  tight,  high  density  implosions  on  axis. 
High  density  is  important  if  efficient  x-ray  production  is  to  be  achieved.  One  of  the  impediments 
to  driving  symmetrical  implosions  is  the  finite  number  of  return  current  posts  that  have  been  used 
to  date.’  They  automatically  destroy  cylindrical  symmetry  and  may  support  the  generation  of 
plasma  turbulence,  which,  in  turn,  softens  die  pinch  and  makes  x-ray  production  more  difficult. 
Current  and  load  asymmetries  will  also  lead  to  2-D  plasma  flow  effects  and  to  Rayleigh-Taylor 
instabilities. 

Whether  or  not  and  by  how  much  load  instabilities  and  asymmetries  degrade  radiation  yield  is 
an  open  question.  2-D  fluid  calculations  are  useful  for  studying  2-D  flow  effects  in  large  aspect  ratio 
implosions  and  for  determining  condititms  for  increasing  imiflosion  stability  and  on-axis  symmetry 
widi  structured  large  diameter  loads.  2-D  fluid  calculations  might  also  be  of  use  in  initiating 
investigations  of  plasma  turbulence,  a  process  not  present  in  1-D  fluid  calculations.  However,  3-D 
calculations  are  needed  if  one  is  to  accurately  track  the  late  time  evolution  of  turimlence.  They 
are  also  needed  to  accurately  investigate  the  nonlinear  evolution  of  Rayleigh-Taylor  instabilities 
and  their  transition  to  turbulent  flow.  If  large  radius  implosions  are  turbulent  (from  finite  numbers 
of  return  currait  paths,  current  driven  vorticity  generation,  or  Rayleigh-Taylor  instabilities),  then 
this  turbulence’s  main  effects  are  already  being  modeled  in  our  1-D  calculations  through  the  use 
of  multipliers  for  viscosity,  heat  conductivity,  and  electrical  resistivity.  An  important  question  for 
large  radius  implosions,  therefore,  is:  will  Rayleigh-Taylor  cause  2-D  average  flow  defects  to  form 
and  how  deleterious  will  tiiey  be  to  x-ray  production?  However,  if  these  effects  can  be  studied  and 
understood,  then  th^  can  be  mitigated. 

20JO  3Q  KsVPhotgns 

There  are  three  basic  processes  by  which  plasmas  radiate;  bound-bound,  free-bound,  and 
free-free.  The  latter  is  also  known  as  bremsstrahlung  emissitm.  Bound-bound  emission  comes  in 
two  forms:  line  emissitm  from  valence  electrons  ot  line  emission  from  inner-shell  transitions.  The 
relative  strength  of  these  processes  in  a  PRS  plasma  is  valence  bound-bound  >  >  firee-bound  >  > 
free-free  ~  inner-shell  line.  One  example  of  free-bound  emission  is  shown  in  Fig.  (14),  which 
contains  a  power  q)ectrum  for  aig<m  tiiat  was  obtained  during  the  course  of  a  1-D  MHD  implosion 
calculatimi.  Notice  that  die  K-shell  lines,  loc^d  at  3  to  4.5  keV  for  argon,  at  on  top  of  two 
free-bound  continue.  Also  notice  tiiat  tiie  continuum  formed  by  the  recombination  of  a  totally 
ionized  atgtm  ion  to  the  hydrogen-like  ground  state  extends  to  10  keV  in  Fig.  (14).  A  dashed  line 
is  drawn  in  tins  figure  to  indicate  that  if  the  argrm  plasma  had  been  hotter  (and  confined),  this 
firee-bound  ccmtinuum  would  extend  to  20  and  30  keV  and  beyond.  Thus,  a  hot,  dense,  confined 
argon  pliama  would  be  an  ideal  source  for  20  to  30  keV  ftiiotons  for  NWET  (on  DECADE?). 

Pgdiatfve  Collapse 

The  K-shell  scaling  laws,  developed  to  date,  make  use  of  1-D  hydrodynamics  calculations 
in  which  the  current  is  terminated  whm  early  slug  model  calculations  wore  tmninated.  These 
cunent-ofif  have  the  utility  of  making  kinetic  energy  the  principle  energy  iiqnit  to 

tiie  plasma  so  tiiat  tiie  calculated  x-ray  conversitm  efficiencies  pertain  only  to  kinetic  energy 
conversion.  In  eaqieriments,  however,  the  current  does  not  turn  off,  and  the  present  scaling  laws 
may  be  too  conservative.  Maintaining  the  current  cm  increases  the  confinement  of  the  pindi,  and 
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it  can  also  lead  to  its  radiative  coUi^se.  Both  of  these  phenomena  could  enhance  the  radiative 
output.  They  may  also  enhance  production  of  x  rays  in  the  20  to  30  keV  energy  range. 

Figs.  (IS)  and  (16)  show  the  relevance  of  current-on  calculations  to  DECADE.  Because  of  its 
large  current,  DECADE  will  have  the  ci^ability  of  accelerating  a  wider  range  of  aluminum  masses 
into  the  efficient  scaling  regime,  depending  on  the  initial  radius,  r^,  of  the  aluminum  array,  than 
existing  simulators  have.  This  point  is  illustrated  by  the  shaded  region  drawn  in  Fig.  (15).  The 
vertical  line,  drawn  at  0.8  mg/cm  as  an  example,  shows  how  DECADE  will  provide  a  capability  to 
mi^  out  the  dependence  of  K-shell  yield  scaling  on  rj  for  fixed  m.  This  capability  is  not  available 
on  curroitly  operational  PRS  simulators.  Fig.  (IS)  also  suggests  that  DECADE  should  provide  a 
capability  of  exploring  radiative  collrpse  behavior  over  a  range  of  masses  from  0.4  to  8  mg/cm. 

The  role  of  kinetic  energy  in  the  implosion  dynamics  will  also  be  much  more  varied  in 
DECADE  than  for  existing  PRS  simulators.  Fig.  (16)  illustrates  this  point.  This  figure  is  to 
Fig.  (IS)  as  Fig.  (3)  is  to  Fig.  (1).  The  numbers  that  label  the  boundaries  in  Figs.  (IS)  and  (16) 
are  in  cme-to-one  correspondence.  Thus,  the  implosions  between  0.2  and  0.6  mg/cm  that  lie  along 
boundary  1  in  Fig.  (IS)  convert  only  a  small  fraction  i'^2  to  4%)  of  the  Marx  bank  energy  into 
kinetic  energy.  A  substantial  portion  of  the  Marx  energy  may  be  available,  therefore,  to  drive 
implosions  that  radiatively  colltpse  on  axis.  DECADE’S  PRS  capability  will  be  essential  for 
carrying  out  these  experiments.  They  could  enhance  simulation  fidelity  capability. 
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Fig.  8  Contours  of  constant  argon  K-shell  3rield  as  a  function  of  the  implosion  velocity  and  mass 
loading  per  unit  length.  Optimal  experimental  results  for  DOUBLE  EAOLE  and  SATURN 
are  plotted. 
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Fig  9  Transmission  line  model  used  for  ACE  IV.  The  pinch  length  is  taken  as  3  cm. 


Current  profiles  for  various  values  of  in  units  of  mg  cm^.  t  =  0  corresponds  to  the 
opening  of  the  plasma  opening  switch. 


(ulo)  “a 


Pig.l  1  Contours  of  constant  ai]gon  K-shell  yield  (dotted)  as  a  function  of  the  initial  puff  radius  and 
mass  loading.  Contours  of  constant  7f  are  shown  as  soUd,  and  constant  in  units  of  mg 
cm^  as  dashed. 
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n.  ONE  -  DIMENSIONAL  MODELING  OF  Z-PINCH  IMPLOSIONS 


One  of  the  major  tasks  of  our  branch  is  to  provide  continuing  data  analysis  and  guidance  for 
z-pinch  experiments  performed  on  DNA  simulators.  A  substantial  part  of  this  analysis  has  been 
facilitated  by  comparing  aq)eiimental  results  witii  predictions  for  the  plasma  <fynamics  that  are 
based  cm  one-dimensional  (l-D)magnetohydro<tyiuunic  (MHD)  calculations.  The  1-D  models  are 
very  important  because  they  present  a  reastmable  tradeoff  between  tiie  two  goals  of  (1)  simplicity 
and  (2)  s(^)histication  -  enou^  to  model  tiie  problems  of  interest.  Specifically,  they  provide 
important  radial  information  (gradients  in  temperature,  vdodty,  density,  radiation,  etc.)  diat  is  not 
present  in  simpler  and  less  costly  zero-dimensitmal  (0-D)  l-cdl  models,  udiich  assume  uniform 
idasma  conditions.  As  an  example,  because  of  tiieir  lack  of  radid  information  0-D  models  are 
unaUe  to  diow,  as  the  1-D  models  are,  tiiat  it  is  generally  tiie  core  and  interior  regimis  of  tiie 
plasma  tiiat  generate  most  of  tiie  K-shdl  radiatitm. 

Some  problems,  such  as  zippering,  nozzle  dynamics,  and  asymmetrical  instability 
development,  requite  multi-dimensional  motkling  tiieir  very  nature.  Unfortunately,  tiie 
models  needed  for  tiiese  problems  either  do  not  exist  or  tiiey  are  so  costly  that  cme  most  usually 
sacrifice  accora^  in  radiation  tran^tt  in  order  to  lower  the  cost  of  induding  extra  dimensions. 
Because  of  tiw  high  cost,  whenever  tiie  plasma  motitm  is  on  average  <nie  dimensitmd  it  is 
advantageous  to  use  a  1-D  modd. 

Our  1-D  models  contain  (1)  a  cmnptehensive  desctiption  of  the  ionization  (tynamics  of 
z-pinch  plasmas,  (2)  an  accurate  description  of  the  tranqpott  of  radiation  within  tiiese  (Really 
thick  plasmas,  and  (3)  the  sdf-ctmsistent  coupling  of  these  dynamics  to  each  other  and  to  the 
fluid  dynamics.  Even  with  tiiese  extensive  equabilities  early  cdculations^’^  predicted  roughly 
two  or  more  orders  of  magnitude  larger  ion  densities  <m  axis  and  coneqxmdingly  much  shorter 
radiation  pulse  widths  than  were  seen  in  experiments.  Because  of  the  lack  of  agreemoit 

between  cdculations  and  eqietiment  we  have  devdoped  over  the  course  of  tiie  last  two  years  a 
phenomenological  modd  for  plasma  turbulence  tiiat  provides  much  better  agrement  between  the 
esqierimnotd  and  cdculated  state  of  tiie  stagnation  plasma  (see  Refs.  4, 8,  and  9).  A  summary  of 
tins  work  is  given  bdow. 
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Summary  of  Phenomenological  Modeling  of  Turbulence 

Based  on  the  lack  of  agreement  between  calculations  and  experiment  and  on  2-D  MHD  z-pinch 
calculations^°’^^*^^  which  showed  the  possibility  of  large  2-D  flow  effects,  it  seems  likely  that 
smne  of  the  missing  |fliysics  in  the  1-D  description  of  the  z-pinch  implosion  dynamics  is  the  onset 
of  a  3-D  hydromagnetic  tuibulence.^^'^^  In  this  case,  a  1-D  calculation  with  modified  transit 
coefficients  might  more  realistically  describe  the  implosion  dynamics  provided  the  turbulmce  were 
uniform  and  the  implosion  was,  on  average,  1-D. 

To  test  this  hypothesis,  we  employed  a  plrenomenological  qrproach  in  order  to  incorporate 
turbulence  effects  into  our  1-D  MHD  calculations.  Multipliers  were  introduced  for  the  plasma 
(artificial)  viscosity,  the  heat  conductivity,  and  the  dectrical  resistivity.  It  was  then  possible  to  make 
I  an  qtpropriate  choice  for  these  multipliers  so  that  die  calculations  produced  qtproximately  the  same 
plasma  conditions  at  stagiuuion  as  a  Physics  International  aluminum  wire  array  experiment.  In 
particular,  comparisons  between  experiment  and  calculations  were  made  for  the  average  stagnation 
ion  densities  and  electron  temperatures  as  well  as  for  the  size  of  the  K-shell  emission  region  and 
the  mass  fractitm  of  the  plasma  emitting  in  the  K-shell. 

We  found  that  it  is  necessary  to  have  large  multipliers  on  viscosity  and  thermal  conduction  in 
order  to  achieve  reasonable  agreement  with  tire  experimental  stagruuion  conditions  of  a  variety  of 
bodi  argon  gas  puff  and  aluminum  wire  implosions.  Large  multipliers  on  viscosity  are  effective 
in  lowering  the  stagnation  ion  densities  and  large  multipliers  on  thermal  conduction  are  effective 
in  convecting  energy  out  of  the  hot  core,  thereby  heating  more  plasma  to  temperatures  that 
are  conducive  for  K-shell  onission  while  simultaneously  lowering  the  core  temperature.  It  is 
difficult  to  determine  a  unique  set  of  multipliers  that  gives  the  best  rq>resentation  for  all  of  the 
experimental  parameters,  but  it  was  found  that  the  calculated  stagnation  conditions  do  ippear  to 
be  smnewhat  insensitive  to  the  choice  of  multipliers  as  long  as  large  multipliers  are  selected.  In 
particular,  the  choice:  40  for  viscosity,  30  for  heat  conduction,  and  20  for  resistivity,  provided 
reasonable  agreement  wifli  all  of  these  experiments,  and  substantially  better  agreemem  than  had 
been  achieved  in  the  non-turbulent  calculations  (see  Figs.  1-3).  This  agreement  would  appear 
to  provide  sipporting  evidmce  for  the  presoice  of  plasma  turbulence  in  z-pinch  expoiments. 
This  phenomenological  approach  to  the  modeling  of  turbulence  in  z-pinch  implosions  based  on 
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benchmaiks  to  experimental  plasmas  may  provide  useful  guidance  for  future  work  that  involves 
tlw  modeling  of  turbulence  in  these  implosions  from  first  principles. 

If  z-pinch  radiation  sources  are  scaled  to  higher  current  machines,  then  it  will  be  important  to 
know  how  turbulence  jdiaiomena  in  these  plasmas  will  scale  as  well.  The  simplest  assumption  to 
make  is  that  it  will  be  as  strong  in  higher  current  machines  as  it  is  in  present-day  machines.  For  this 
reason,  we  recalculated  die  yield  scaling  with  mass  for  both  aluminum  and  argon  plasmas  using  the 
same  transport  coefficient  oihancements  as  had  been  determined  in  the  analysis  of  the  PI  aluminum 
and  argon  experiments.  Softening  the  implosions  has  the  consequence  of  requiring  substantially 
more  mass  (factor  of  6)  to  achieve  efficient  K>shell  yield  scaling  (for  t;  =  4)  (see  Figs.  4  and  S), 
since  turbulence  lowers  the  stagnation  ion  dcmsity  for  a  given  mass  loading.  Consequendy,  the 
ability  of  the  plasma  to  radiate  is  diminished  for  a  given  mass,  and  the  radiative  cqiacity  can  be 
increased  only  by  imploding  more  mass.  This  factor  of  6  has  been  of  major  importance  for  making 
future  predictions  about  the  K-shell  yield  petfonnance  of  DECADE  and  JUPITER  simulators. 

Future  Directions  of  1>D  Z-Pinch  Modeling 

From  our  work  on  {dienomenological  modeling,  which  includes  making  numerous 
(xmiparisons  between  experimental  and  calculated  results,  we  know  that  there  is  still  a  great 
deal  to  learn  and  understand  about  z-piach  behavior.  It  will  be  necessary  to  obtain  this  knowledge 
and  understanding  in  order  to  design  the  plasma  radiation  source  (PRS)  loads  tiiat  DNA  requires 
of  its  future  simulators.  Below  we  discuss  a  few  of  the  areas  of  z-pinch  research  that  we 
are  investigating  with  our  1-D  modelling  ciq)abilities.  Note,  there  are  many  multi-dimensional 
problems,  such  as  optimal  nozzle  design  and  minimizing  instability  growth  ,  that  are  described  in 
otlwr  sections  of  this  report. 

(a)  post  implosion  modeling  and  short  circuiting 

For  implosions  that  occur  late  in  the  current  pulse  and  for  those  machines  which  cannot  put 
much  current  into  the  load  at  stagnation  because  of  large  induction,  late  time  currait  flow  and 
confinement  may  not  be  so  important  for  determining  K-shell  emission.  However,  for  several  of 
tiie  Double  Eagle  aigtm,  krypton,  and  aluminum  e3q)eriments  that  produced  substantially  more 
K-shell  emissinn  than  predicted  by  current-off  turbulence  calculations,  it  is  likely  that  an  accurate 
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model  of  the  late  time  current  effects  is  needed.  The  straightforward  tq^proach  of  just  allowing 
the  current  to  flow  during  stagnation  produces  results  which  are  not  generally  in  accord  with 
experiment.  For  example,  calculations  made  for  low  mass  loads  such  as  the  42  tig/ cm  I%ysics 
International  (PI)  e:q)erimental  argon  load  show  that  keeping  the  current  on  leads  to  K-shell  yields 
mote  in  accord  with  experimoit  but  it  also  leads  to  an  elastic  bouncing  behavior  that  causes  the 
K-sheU  radiatitm  to  be  emitted  over  100  ns  (see  Fig.  6).  Neither  the  bounce  behavior  or  the 
l(Hig  emission  times  are  characteristic  of  the  experiment.  One  can  q)eculate  that  part  of  long 
emission  time  behavior  could  be  due  inherendy  to  die  1-D  nature  of  the  model,  which  only  allows 
uniform  radial  modtm.  If  the  theoretical  plasma  is  allowed  to  move  axially  after  sta^iation,  then 
part  of  the  stored  internal  energy  will  not  be  available  to  do  work  against  the  magnetic  field.  The 
result  will  be  a  smaller  radial  expansion  and  a  shorter  time  between  radial  contractions. 

As  another  example,  current-on  calculations  made  for  the  more  massive  PI  aluminum  loads 
(474, 277,  and  164  yg/cm)  produced  radiative  coll^>ses  for  which  there  was  more  energy  coupled 
into  the  load  and  radiated  in  the  K-sheU  than  was  observed  in  the  experiments.  For  this  reason  we 
will  omtinue  to  investigate  the  possibility  of  alternative  current  paths  and  shorting  of  the  diode 
which  could  reduce  the  energy  coupling  to  tte  load  and  magnetic  field  confinement  of  the  plasma. 
Some  of  this  wodt  as  it  relates  to  Gamble  n  and  Double  Eagle  was  presented  in  the  1990  DNA 
annual  report.  More  recently,  a  short  circuit  analysis  of  the  Saturn  aluminum  experiments  showed 
that  the  total  radiative  yield  behavior  that  was  exhibited  in  the  e:q[)eriments  had  the  same  trend  as 
the  energy  cotq)led  into  a  shorted  diode.  The  following  example  is  for  a  short  that  occurs  when  the 
voltage  across  the  diode  is  3MV. 


Mass 

/ig/cm 

Anode 

radius 

(cm) 

Total  Exp. 
Radiation 
(kJ) 

Coupled  Energy 
3MV  short 

OJ) 

K-sheU 

Yield 

Exp.-(kJ) 

K-sheU 

Yield 

Theory -(kJ) 

120 

2.4 

172 

235 

8 

14 

160 

2.4 

170 

258 

12 

20 

185 

2.4 

196 

278 

22 

24 

330 

2.4 

331 

346 

63 

79 

470 

2.4 

364 

383 

61 

120 

740 

2.4 

338 

433 

30 

76 
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While  this  work  does  not  suj^ly  1 00%  conclusive  evidence  that  shorts  are  playing  a  predictable 
role  in  the  tfynamics  of  the  implosicms,  it  does  dem<»strate  that  shorts  deserve  more  experimeittal 
uid  theoretical  ccmsideration. 

(b)  generating  >  SO  keV  photons 

Our  estimates,  which  are  based  <m  qjjdying  cuiient  K-shell  yield  scaling  results  to  xenon 
(ZsS4)  plasmas,  show  tiiat  at  best  tots  of  Idlojmiles  of  30  keV  phoUHis  mi^  be  obtained  using 
Jiqater  class  machines  (see  .^ipendix).  To  address  the  problem  of  attaining  significant  yidds 
(hundreds  of  Idlojoules)  of  >  30  keV  idx>ti)n8  using  future  simulators,  we  began  an  investigation 
into  examining  the  efficiency  ofproducing  high  energy  recombination  |Aotong.  This  scheme  entails 
using  a  z-pinch  to  heat  the  core  of  a  moderate  Z  material,  Ni  for  example,  to  very  hot  temperatures 
and  then  rely  on  radiative  reoomlanation  to  the  hdiumlike  and  lydrogenlike  states  to  produce  the 
high  energy  photons.  Thermal  bremsstrahlung  is  also  present  but  its  coefficient  is  substantially 
smaller  tiian  that  fig  recombination  radiation.  Since  the  fipee-bourxl  (and  Bremsstrahlung)  emissitm 
is  proportional  to  Z^Ni^exp{=j^),  tiiere  ritould  be  proportionately  more  i^iotons  coming  out  at 
hi^  erietgy  as  tiie  core  temperature  of  tiie  jdasma  is  raised.-  We  do  not  expect  to  attain  efficient 
yield  scaling  of  30  keV  jrfiotons  using  tiiis  scheme;  however,  we  do  fed  diat  we  can  do  better  tiian 
would  be  attained  by  im|tioding  xenon  on  Jiqnter.  The  major  reason  bdiind  tiiis  qieculatitm  is  tiiat 
the  stagnation  densi^  (i\rt)is  much  larger  (udiich  implies  tiie  radiation  effidency  can  be  larger) 
when  imploding  large  mass  moderate  z  plasmas  as  opposed  to  low  mass  hi^  z  plasmas. 

The  study  we  are  initiating  involves  fixing  tiie  mass  of  a  riickd  plasma  at  2000  fig/cm  and 
fixing  tiie  initial  radins  at  2.0  cm.  The  load  is  tiien  imitioded  using  a  linearly  ramped  currem  profile 
tiiat  is  terminated  when  the  outer  radius  of  tiie  pinch  reaches  0.2  cm.  Under  these  conditions  we 
found  tiiat  tiie  >  30  keV  radiation  is  a  maximum  for  a  currem  profile  which  produces  a  peak 
currem  of  51  MA  in  50.4  ns.  Ihe  work  done  on  this  plasma  is  =  4  MJ/cm,  whidi  is  wdl  within 
tiie  realm  of  a  Jupiter  class  machine.  Fig.  7  shows  what  tiie  predicted  yidds  are  above  each  energy 
levd  for  the  51  MA  calculation  and  in  particular  it  shows  tiiat  27  kJ/cm  of  >  30  keV  photons  are 
being  produced.  Since  we  have  only  sampled  a  small  portion  of  the  available  parameter  space  in 
terma  of  mass  loading  and  energy  coupling  to  the  load  tiiete  is  good  reason  to  be  optimistic  that 
more  than  27  kJ/cm  of  >  30  keV  radiatimi  can  be  goierated  Jiqnter  dass  machines.  We  are 
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also  investigating  DECADE’S  ciq>ability  of  producing  >  20  keV  photons. 

(c)  diagnosing  core  physics 

Because  our  1-D  model  uses  cell-to-cell  coupled  probability  of  escape  techniques  to  transport 
radiati(m^‘*^*'^^  there  is  an  inherent  uncertainty  in  inferring  the  final  emitting  zoim  of  an  escq>ing 
photcm  because  of  the  averaging  over  the  line  profile.  Thus,  q>atial  comparisons  between  the  sizes 
of  die  experimental  emission  regions,  which  are  based  mi  x-ray  pinhole  pictures,  and  our  calculated 
sizes  are  smnewhat  problematical.  This  has  not  been  a  problem  for  analyzing  die  PI  eiqietiments 
because  it  was  usually  only  the  inner  regions  of  the  plasmas  diat  were  hot  enmi^  to  itmize  into 
die  K-shell  udiile  die  outer  regions  remained  too  cold  and  underimiized  to  be  efifecdve  reabsorbers 
of  K-shdl  radiadon.  Therefore,  even  thou^  the  exact  zone  from  which  a  K-diell  {diotmi  escqied 
was  unknown,  the  size  of  die  sonroe  region  of  plasma  from  v^ch  it  esc^ied,  in  diese  calculations, 
was  faidy  well  defined.  However,  vhen  we  began  an  analysis  of  die  Saturn  aluminum  data  we 
found  evidence  diat  dus  is  not  die  case.  For  example  vdien  we  die  qpectioscopic  techniques 
of  Coulter  etal.  to  the  Saturn  data  we  find  that  tlK  average  temperature  of  the  K-shell  emissicHi 
region  of  die  best  K-diell  emitting  load  (330  pg/em)  is  600  -  900  eV  widi  approximately  30  % 
of  the  mass  patdcqiating  in  producing  K-shell  emissimi.  The  problem  is  dut  die  600  -  900  eV 
temperatures  seem  to  be  mudi  lower  dian  expected  because  we  know  dua  similar  mass  loadings 
for  argrni  have  produced  consideraUe  K-shell  emission  on  Saturn,  vdiich  would  require  much 
higher  temperatures.  This  is  also  in  accord  widi  our  1-D  analysis  which  shows  the  average  K-shell 
emission  regirm  tmnperature  to  be  ipproximately  1600  eV.  It  is  also  in  accord  widi  an  earlier  study 
in  which  wc  if^ed  our  qiectroscopic  technique  to  a  1-D  aluminum  calculation  and  found 
that  the  ^lectroscopic  averaging  process  predicted  smaller  temperatures.  At  diis  point  we  are 
investigating  the  possibility  that  there  is  substantial  absorption  and  reemission  of  K-shell  radiation 
taking  place  in  cooler  regions  of  the  Saturn  plasmas  that  are  affecting  the  Lyman  and  He-alpha  line 
ratios  and  consequoidy  our  qiectroscqpic  analysis. 

In  the  past  we  have  always  distributed  some  mass  of  the  plasma  all  die  way  to  die  axis. 
This  procedure  has  provided  numerical  stability  by  providing  back  pressure  early  in  the  stagnation 
diat  softens  die  implorinn.  This  procedure  may  be  adequate  for  gas  puff  simulations  because 
of  die  divergence  of  the  gas  from  die  nozzle  towards  the  axis,  however,  it  is  not  clear  that  wire 
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array  implosicms  behave  in  this  same  manor.  Therefore  we  have  begun  an  investigation  into  the 
differences  in  stagnation  dynamics  betwem  a  shell  implosion,  which  may  be  more  indicative  of 
a  wire  array  experiment,  and  a  snow  plow  imf^ositm,  which  is  more  characteristic  of  a  gas  puff 
calculation. 

Appendix:  Simple  K>Shell  Yield  Scaling 

After  ctmversatiions  widi  members  of  die  PRS  community  it  became  iqipaient  diat  several 
members  were  having  difficulty  in  qiplying  our  PRS  K-shell  yield  scaling  results  to  DECADE 
and  JUPITER  class  machines,  ^ipaiendy,  d»  community  is  interested  in  a  singular  choice  of 
initial  load  conditions  diat  will  t^itimize  K-shdl  emission,  whereas,  our  results  predict  that  there 
is  a  l»oad  range  cf  mass  loads  and  initial  radii  diat  will  produce  neatly  die  same  t^itimal  K-diell 
yield  from  a  given  dement  and  machine.  In  order  to  facilitate  usage  of  our  woik  by  the  PRS 
community  we  dedded  to  pnqiose  our  own  simple  K-slidl  yidd  scaling  rdation.  This  rdation 
makes  a  singular  judidous  choice  of  initial  load  amdititms  diat  is  founded  upon  the  prindples  of 
our  eariier  1-D  scaling  woik.  Based  cm  die  kinetic  energy  coupling  E  firmn  the  machine  to  the 
load  dus  scaling  relation  predicts  die  qitimal  mass  load  M«pt  and  the  expected  K-shell  emission 
as  a  function  of  atomic  number  Z  of  die  load.  Hie  intid  radius  can  be  found  using  a  simple  slug 
modd  or,  if  it  is  i^ipticable,  a  sinusoidd  modd  can  also  be  used.  The  scaling  rdatim  is  derived 
cdcularing  die  largest  vahie  of  defined  as  diat  can  be  achieved,  for  a  given  machine 
kinetic  energy  coupling  to  the  load  E  (j/cm),  vdiile  still  imploding  enough  mass  to  insure  effici<mt 
K-diell  yidd  scaling.  That  is  i7m«c  satisfies  the  equatitm 

n  Emin  Mbp{v)I^  =  ^  (1) 

where  A  is  the  atomic  mass  (g),  Emin  is  the  mininmim  energy  needed  to  ionize  into  die  K-shell 
(j/ion)>  and  Mbp  is  die  break-point  mass  (g/cm).  All  of  diese  quantities  are  defined  in  Ref.  2. 
Note  Mbp  has  been  idieiiomenologically  determined  to  be  a  factor  of  6  greater  than  the  vdue 
quotedinRef.  2.'  A  useful  ai^iroximation  for  the  kinetic  energy  coupling  E  is  given  by 

£  =  1607/*  (j/cm)  (2) 

vrhae  I  is  die  peak  load  current  (MA).  Equaticm  (1)  can  be  solved  numerically  for  rjmaz  or  it  is 
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given  i^TOximately  by 


I/m..  S  (6.3  /or  1/^.  <7.  (3) 

The  K-shell  yield  scaling  Yk  (j/cm)  is  given  by: 

Yk  =  0.3E  (j/em)  for  17^,  >  2.0  (4) 

y*  =  0.3£(-i/^.  +  Arimaz  -  3)  1-25  <  »/m..  <  2.0 

The  (^)timal  mass  loading  -  MBp(f/m..)  and  it  is  iqiproximately  given 

s  10-T  {g/em)  for  1.25  <  i/«.,  <  7  (5) 

For  t/m..  <  1.25  the  following  relations  zpply 

=  1.3  10-®  {g/em)  (6) 

n  =  1.3  10'®JlV^'*(i/cm)  (7) 

Eqiiatieas'(6)  and  (7)  draw  that  a  Jupiter  machine  coiqding  S  MJ/cm  of  kinetic  energy  to  a 
load  can  expect  to  adiieve  17  kJ/cm  of  xenon  K-shell  emission.  Note,  there  is  still  a  great  deal  of 
uncertainty  in  die  yidd  scaling  at  low  17. 
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<Ion  Densjty>  (10^®  cm" 


Figure  1.  Laminar,  t  jbulent,  and  experimental  stagnation  ion  densities  as  a  flmction  of 
almniniitn  load  masS. 
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<Electron  Temperature>  (eV) 


Figure  2.  T  -aminiir,  tuibulcot,  and  experimental  stagnation  electron  temperatures  as  a  ftmction 
of  aluminum  load  mass. 
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Mass  (;ig/cm) 


Figure  3.  Laminar,  turbulent,  and  caqMiimental  mass  fractions  of  die  plasma  emitting  in  the 
K-diell  as  a  function  of  aluminum  load  mass.  The  multipliers  on  the  transport  coefhdents  for 
die  turbulent  calculations  are:  30,  40,  and  20  for  beat  conductivity,  viscosity,  and  reasdvity, 
respectively. 
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Figure  4.  Non-tuibulent  aluminum  K-shell  yield  scaling  with  mass  for  17  =  4. 
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Figure  5.  *nut>uleot  aluminum  K-shell  yield  scaling  with  mass  for  17  s  4.  The  multipliers 
on  the  transport  coefiGicients  are:  30, 40,  and  20  for  heat  conductivity,  viscosity,  and  resistivity, 
respectively. 
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Figure  6.  Hie  calculated  time  evolution  of  K-sbdl  emission  power  firmn  a  42  yglem  load  is 
diown.  b  tile  calculation,  the  cunent  and  enhanced  transport  are  kept  on  throughout  several 
stagnation  and  rebound  t^des  of  tiiepindL 
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Figure  7.  Tune  integrated  emisaon  above  each  {dioton  energy  levd  for  a  2000  fig /cm  nick<d 
load. 


ni.  PRS  Iioad  Simulation 


Throughout  this  past  year,  several  topics  in  the  area  of  2-D  modeling  were 
addressed:  (1)  the  viability  of  PRS  loads  in  large  current  machines  (i.e.,  on  the 
order  of  30  MA  and  larger),  (2)  Support  of  ongoing  Gas  Puff  Experiments,  and  (3) 
code  improvement  and  modification. 

(1)  Large  Current  Machines 

During  the  past  year,  a  large  number  of  simulations  were  made  which  aimed 
at  evaluating  PRS  loads  on  a  hypothetical  large  current  driver.  We  concentrated 
on  80  MA,  60  MA,  and  40  MA  peak  current  in  the  load  during  the  runin  from  an 
initial  radius  of  8  cm  (a  3cm  diode  length  was  assumed).  These  were  Krjrpton  gas 
puff  loads.  The  driving  current  in  these  simulations  used  the  circuit  model  with  a 
70  nsec  switch  time  and  a  12  nh  inductor  on  the  inductive  store  side  of  the  circuit. 
An  additional  assumption  was  made  that  there  is  a  total  of  2.5  nh  inductance  on 
the  load  side  of  the  switch  which  is  coupled  to  an  additional  2.5  nh  short  circuit 
load  inductance.  Naturally,  this  last  source  of  inductance  increases  as  the  implosion 
proceeds  and  the  radius  decreases.  The  initial  energy  in  the  storage  inductor  was 
adjusted  until  the  desired  load  current  during  the  runin  was  achieved.  The  results 
from  our  2-D  computer  code,  PRISM,  operating  in  a  1-D  mode  show  peak  K- 
shell  yields  of  approximately  860  kJ/cm,  400  kJ/cm,  and  100  kJ/cm  for  80  MA, 
60  MA,  and  40  MA  peak  load  currents.  The  yields  were  obtained  for  implosions 
occurring  at  140  nsec  and  masses  which  were  3.3  mg/cm,  2.3  mg/cm,  and  1.1 
mg/cm,  respectively.  A  caveat  for  these  results  is  the  fact  that  a  radiation  lookup 
table  routine  with  no  transport  was  used.  This  can  overestimate  the  radiative  yield 
when  larger  masses  are  used.  For  instance,  1-D  results  at  higher  mass  levels  show 
yields  which  decreased  unrealistically  as  the  mass  was  increased. 

Our  subsequent  work  on  2-D  modeling  have  shown  that  loads  which  start  at 
this  large  a  radius  are  highly  unstable  (see  discussion  below).  The  requirement  for 
uTTiiJW  radii  greatly  reduces  the  kinetic  energy  which  can  be  converted  to  radiation. 

One  of  the  important  issues  for  large  radii  loads  is  the  relevance  of  the  Rayleigh 
Taylor  (R-T)  instability.  We  have  modeled  this  instability  in  two-dimensions  using 


41 


parameters  relevant  to  PRS  loads  on  future  pulse-power  generators.  We  examined 
implosions  from  a  10  cm  radius,  with  a  3  mg/cm  Krypton  mass  loading,  and  a  50 
MA  current  in  the  load.  In  these  calculations  a  moving  mesh  was  used  with  the 
highest  resolution  concentrated  in  the  vicinity  of  the  peak  density  of  the  compressed 
shell.  The  mesh  in  these  simulations  was  continually  adjusted  to  remain  at  10  cm  by 
following  the  peak  in  the  shell  density  in  time.  This  simplification  can  be  justified 
because  the  aspect  ratio  (radius  of  the  shell/ width  of  the  shell  1)  is  large.  In 
addition,  this  procedure  allowed  us  to  start  the  2-D  instability  runs  with  steady  state 
profiles  obtained  horn  initial  1-D  runs.  In  this  fashion,  compressive  transients  which 
could  spoil  the  desired  mode  structure  are  avoided.  These  results  have  shown  that 
the  growth  rate  is  on  the  order  of  1-3  x  10^  sec~^  up  to  10*  sec“^  for  a  wavelength 
roughly  equivalent  to  the  shell  thickness. 

Next,  a  series  of  simulations  were  made  which  started  with  initial  1-D  solutions 
but  in  which  the  compressed  shell  was  allowed  to  move  through  the  mesh,  i.e. 
implode.  These  annular  gas  puff  shells  exhibited  Rayleigh- Taylor  instabilities  when 
accelerated  from  initial  radii  of  approximately  10  cm.  The  instabilities  were  found 
to  have  a  growth  rate  as  high  as  1  -  2  x  10*  sec~^  which  amounts  to  an  e-folding  time 
of  5  •  10  nsec.  This  is  in  agreement  with  the  fixed  location  results  discussed  above. 
The  somewhat  higher  growth  rates  come  about  because  during  the  implosion,  the 
acceleration  increases  as  the  pinch  moves  inward.  In  simulations  of  a  gas  annulus 
imploding  from  8  cm,  growth  rates  were  measured  as  high  as  1.3  x  10*  sec~^  for  short 
wavelength  perturbations  (50  MA  peak  current).  Therefore,  only  35  nsec  is  required 
for  a  1  percent  perturbation  to  grow  exponentially  to  100  percent.  Obviously,  this 
is  well  before  peak  radiation  wotild  occur  at  majdmum  compression  on-axis.  This 
disrupts  the  imploding  gas  well  before  the  final  collapse  stage.  In  this  particular 
case,  the  pinch  had  moved  from  8  cm  down  to  about  5-6  cm  before  disruption. 

We  farther  examined  the  stability  of  implosions  with  radii  from  4  cm  up  to 
10  cm.  These  implosions  are  driven  by  a  simple  circuit  (described  above)  which 
consists  of  an  inductor  with  roughly  120  MA  switched  into  the  load  with  a  70  nsec 
switch.  This  results  in  roughly  a  50  -  60  MA  driving  current  during  the  implosion. 
However,  it  is  important  to  note  that  scaling  laws  can  be  applied  so  that  the  resiilts 
should  be  siinilar  for  lower  currents  as  well.  For  instance,  changing  the  mass-radius 
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product,  MR’,  would  result  in  a  lower  peak  current,  /,  and  shorter  implosion  time 
for  the  same  circuit  but  /  would  remain  the  same.  Parameters  within  the  circuit  can 
be  altered  to  adjust  /.  A  simple  calculation  shows  why  such  a  configuration  is  very 
susceptible  to  the  classical  Rayleigh-Taylor  instability  with  growth  rate  7  =  y/leg, 
where  k  =  2ir/X  and  g  is  gravity  or  acceleration.  If  one  assumes  an  approximately 
uniform  acceleration,  a,  the  distance  traveled  by  the  shell  is  pven  as  d  =  where 
t  is  time.  The  amplification  of  the  initial  perturbation  is  then  with  =  ^2hd. 
Note  that  this  term  is  independent  of  both  time  and  acceleration.  It  is  well  known 
that  for  shells  undergoing  acceleration,  the  most  dangerous  growth  occurs  when 
kA  ~  1  where  A  is  the  shell  thickness  which  itself  depends  on  the  shell  mass, 
acceleration,  etc.  For  the  moment,  however,  if  we  take  k=27r/0.2  cm  and  assume 
that  the  shell  is  completely  disrupted  when  a  perturbation  is  amplified  by  a  factor 
of  10^,  the  shell  will  travel  only  about  1.3  cm  before  disruption  occurs.  This  is  an 
idealised  argument  which  assumes  dassical  planar  geometry  with  a  sharp  boundary. 
However,  it  is  realistic  enough  to  shed  light  on  the  difficulty  in  driving  thin  shells  in 
from  large  distances  as  is  being  proposed  for  future  PRS  machines  (i.e.,  DECADE 
and  JUPITER).  This  is  also  consistent  with  the  results  dted  above. 

Simulated  2-D  implosions  using  the  drcuit  outlined  above  with  two  different 
fill  conditions  are  presented  in  Figures  1  and  2.  These  simulations  were  for  initial 
plasma  configurations  which  are  shells  (i.e.  A/R  >>  1)  and  uniform  fills. 

Figure  1  shows  the  time  history  of  the  fraction  of  mass  in  zones  with  greater 
than  1.25  the  average  mass,  the  fourier  component  across  the  colunm  of  the  mass 
difference  for  4  nodes,  and  the  fourier  component  for  2  nodes.  The  mass  of  the  shell 
was  3.33  mg/cm  (Krypton)  in  an  initial  shell  configuration  starting  at  a  radius  of 
6  cm.  The  initial  perturbation  was  introduced  with  a  slight  kinking  of  the  colunm 
at  four  evenly  spaced  nodes.  In  Fig.  1,  it  can  be  seen  that  there  is  linear  growth 
for  the  four  node  mode  followed  by  a  phase  of  more  rapid  growth  during  which  the 
mass  accumulates  rapidly  in  the  Rayldgh-Ibylor  spikes. 

This  accumulation  is  evidenced  by  the  curve  which  shows  that  the  fraction  of 
mass  in  zones  with  greater  than  1.25  of  the  average  mass  grows  rapidly  to  1.  At  the 
end  of  the  nm,  the  disrupted  shell  had  moved  to  a  radius  of  about  4  cm  (i.e.,  over 
a  distance  of  2  cm).  This  is  roughly  the  disruption  distance  pven  by  the  simple 
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Argument  outlined  above  and  is  not  encouraging  for  a  thin  shell  design  which  must 
implode  over  a  large  distance.  Figure  2  shows  the  time  history  (the  same  quantities 
are  plotted  as  in  Figure  1)  but  for  a  uniform  fill  implosion  with  the  same  mass 
and  circuit  as  in  the  shell  implosion.  In  this  case  it  can  be  seen  that  the  initial 
mode  grows  linearly  but  is  then  suppressed  while  the  2  node  mode  grows  and  also 
is  subsequently  suppressed. 


Figure  1  Fraction  of  mass  in  zones  with  greater  than  1.25  the 
average  mass  and,  fourier  components  of  the  mass  variation  for  2 
and  4  nodes  for  a  shell  implosion. 

The  result  is  that  the  mass  accumulation  growth  is  delayed  when  compared  to  the 
shell  implosion.  The  final  state  at  about  112  nsec  consists  of  a  single  node  (i.e.,  1/2 
wavelength)  in  the  computational  box  with  the  dense  portion  of  the  plasma  well 
inside  1  cm.  This  type  of  implosion  offers  hope  that  plasma  may  reach  the  axis  over 
a  long  runin  distance  without  being  totally  destroyed  as  happens  with  a  thin  shell. 
Recently,  Goldberg  and  Velikovich  (Phys.  Fluids  B  5,1164(1993))  have  proposed 
that  an  implosion  of  a  uniform  fill  configuration  can  lead  to  some  suppression  of 
the  Rayleigh-Taylor  instability  through  accretion  of  mass  onto  the  imploding  shell. 
This  work  seems  to  support  the  results  of  our  simulations.  The  penalty  for  such 
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designs,  however,  is  that  in  uniform  fill  loads  a  large  portion  of  the  mass  may  not 
reach  the  required  kinetic  energy  to  significantly  radiate  in  the  K-shell. 

The  task  in  the  upcoming  year  will  be  to  optimize  designs  for  large  radii  PRS 
loads  which  produce  the  necessary  bulk  kinetic  energy  while  maintaining  some  in¬ 
tegrity  during  the  implosion.  This  investigation  could  have  significant  impact  on 
the  design  of  PRJS  sources  for  future  pulse-power  simulators. 

Our  conclusions  to  date  on  the  stability  of  imploding  loads  to  the  Rayleigh 
Taylor  instability  can  be  summarized  as  follows:  Imploding  shells  appear  to  break 
apart  (over  several  cm  of  runin  length)  for  the  cases  we  have  considered:  10  cm 
down  to  about  5  cm.  Note  that  these  are  radius  values,  not  diameters.  A  simple 
argument  based  on  the  acceleration  determined  from  d  =  1  /2a<^  and  a  growth  rate 
given  as  y/ka,  where  k  is  the  inverse  wavenumber,  gives  a  distance,  d,  of  roughly 
2  cm  before  significant  breakup  of  thin  shell  would  occur.  Our  simulations  show 
that  somewhere  around  a  3  cm  radius,  implosions  of  Krypton  gas  shells  would  seem 
to  be  viable.  This  analysis  is  based  on  the  series  of  runs  discussed  above  and  the 
observed  distances  over  which  disruption  occurred. 

Uniform  fill  implosions  present  a  different  story.  We  have  been  able  to  bring 
the  load  in  from  5-6  cm  and  have  observed  an  interesting  phenomena.  That  is,  as 
the  implosion  proceeds,  the  snowplow  picks  up  extra  mass  from  the  fill  gas  which 
broadens  the  shocked  snowplow  repon  and  suppresses  the  initial  wavelength.  The 
result  is  that  eventually  only  the  longest  wavelength  available  in  the  computational 
box  (1/2  wavelength  from  z=0  to  z=zmax)  is  seen.  For  numerical  resolution  reasons, 
the  length  along  z  cannot  bet  set  to  the  size  of  the  machine  because  the  wavelengths 
are  related  to  the  width  of  the  snowplow.  This  is  much  narrower  than  the  diode 
length.  Therefore,  it  should  be  noted  that  this  result  may  be  somewhat  dependent 
on  how  large  we  set  the  size  of  the  numerical  box.  Uniform  fill  implosions  appear  to 
be  more  stable  than  thin  shell  implosions.  On  the  negative  side,  significant  heating 
is  seen  in  these  implosions  ahead  of  the  snowplow  which  leads  to  hollow  implosions. 
Currently,  our  effort  is  being  directed  toward  assessing  the  role  of  radiation  in  these 
implosions.  Analysis  of  the  uniform  fill  runs  show  that  the  implosion  velocities  reach 
5x10^  cm/sec  after  about  65  nsec  of  a  100  nsec  implosion.  Near  final  implosion, 
it  appears  that  veloaties  as  high  as  8x10^  cm/sec  can  be  achieved.  This  should 


provide  enough  energy  to  produce  Krypton  K-shell  radiation. 


Figure  2  Fraction  of  mass  in  zones  with  greater  than  1.25  the 
average  mass  and,  fourier  components  of  the  mass  variation  for  2 
and  4  nodes  for  a  uniform  implosion. 

(2)  Support  of  Ongoing  Gas  Puff  Experiments 

Several  computer  simulations  were  made  in  support  of  planned  and  ongoing 
experiments  on  the  Double-Eagle  pulse  power  generator  at  Physics  International 
(PI).  These  experiments  are  designed  to  test  new  nozzles  with  smaller  exit  aper¬ 
tures.  The  narrowest  aperture  used  to  date  has  been  4  mm.  In  the  new  designs,  this 
will  be  reduced  to  2  mm.  One  of  the  important  areas  in  which  the  NRL  radiation- 
hydrodynamics  branch  can  contribute  in  the  future  is  to  provide  simulation  studies 
of  the  gasdynamics  of  innovative  nozzle  designs.  This  is  more  important  for  up¬ 
coming  machines  where  the  tradeoffs  discussed  in  the  previous  section  requiring 
structuring  the  design  with  multiple  puffs,  for  instance.  This  necessitates  a  clearer 
understanding  of  the  nozzle  gasdynamics.  In  the  past,  both  our  simulations  and 
the  PI  experiments  have  shown  improvement  when  the  aperture  was  reduced  from 
8  mm  to  4mm.  However,  simulations  performed  during  this  past  year  have  shown 
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only  small  differences  in  yield  when  a  further  reduction  is  made  to  2  mm.  Table  I 
shows  the  results  for  the  standard  Double  Eagle  parameters  which  have  been  used 
in  the  past.  That  is,  1.25  cm  radius,  50  /ig/cm  Argon  mass,  and  3.3  MA  peak 
current.  A  linear  current  ramp  model  was  used.  Both  nozzles  were  tilted  inwards 
at  a  -10**  tilt. 

Nozzle  Exit  Aperture  Yield 

0.4  cm  6  kJ 

0.2  cm  5.2  kJ 

'Dible  I. 

Ehcamination  of  the  simulation  results  show  that  thinner  shells  are  more  prone 
to  breakup  on  the  way  m  and  appears  to  be  the  cause  of  the  lower  yield.  This  phe¬ 
nomena  is  witniluT  tothe  results  discussed  in  the  previous  section  and  demonstrates 
its  relevance  to  current  experiments. 

The  experimental  schedule  will  now  take  this  into  account.  Fewer  of  the 
new  nozzles  with  2  mtn  apertures  are  being  constructed  than  had  been  originally 
planned.  However,  there  will  be  a  few  shots  made  with  the  2  mm  aperture  nozzles 
in  order  to  compare  simulation  and  experiment.  This  collaboration  demonstrates 
that  the  cooperation  between  computer  simulations  and  experimental  planning  can 
restdt  in  savings  of  time  and  money. 

(3)  Code  Improvement  and  Modification 

In  support  of  our  2-D  work  on  gas  puff  experiments,  the  computer  code  PRISM 
is  being  made  available  to  several  users.  That  is,  it  is  presently  making  the  transition 
from  a  single  user  environment  to  a  multi-user  production  environment.  This  has 
required  that  the  interface  and  post-processing  packages  be  made  more  user  friendly 
and  adaptable.  In  addition  to  its  use  in  simulating  PRS  load  dynamics,  PRISM  will 
be  used  in  studying  the  dynamics  in  other  areas  of  interest  to  DNA  such  as  switch 
plasmas.  The  code  was  modified  so  that  it  could  run  on  the  Los  Alamos  DNA  Cray 
machines  under  UNICOS  as  well  on  large  UNIX  based  workstations. 
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(4)  Papers  presented  on  this  work 


A  paper  on  our  previous  work  dealing  with  issues  in  1-D  and  2-D  PRS  implo¬ 
sions  was  presented  at  the  Third  International  Conference  on  Dense  Z-pinches  at 
Imperial  College,  London.  The  abstract  reads:  ‘^Non-LTE  radiation  hydrodynamic 
numerical  simulations  in  1-D  nad  2-D  are  performed  for  multi-terawatt  driven  argon 
and  krypton  gas  puff  loads.  The  influence  of  enhanced  transport  coefficients  on  the 
plasmas’  implosion  dynamics  and  the  effects  it  has  on  the  radiation  yield  and  spec¬ 
tral  distribution  are  discussed.  Also  presented  are  simulations  for  the  performance 
of  a  krypton  gas  puff  driven  by  a  class  of  future  drivers.” 

Some  of  the  work  discussed  above  was  presented  in  an  invited  talk  at  the  IEEE 
meeting  in  early  June  on  the  physics  of  PRS  loads:  Radiation  Hydrodynamics  of 
Z-Pinch  Plasmas,  J.  Davis,  IEEE  Conference  7-9  June  1993,  Invited  Paper  4B5-6. 
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IV.  DECADE  Power  Flow 


•  Coupling  Issues  for  DECADE 

The  introduction  of  inductive  energy  stores  (IBS)  to  pulse  power  driven  Plasma 
Radiation  SourcesCPRS)  has  made  possible  an  order  of  magnitude  or  so  reduction  in  the 
total  size  of  the  pulser  modules  and  a  potentially  significant  decrease  in  the  risetime 
of  the  output  current  pulses.  NRL  is  involved  in  a  specific  effort  to  understand  power 
flow  to  the  PRS  load  in  an  IBS  machine  in  support  of  the  DNA  DBCADB  simulator. 
As  part  of  that  program  we  have  developed  and  tested  a  transmission  line  model  that 
resolves  the  Plasma  Opening  Switch  (POS),  and  connects  that  model  to  a  simple  "gas 
bag"^^^  PRS  model.This  work  is  intended  to  provide  a  simplified  but  comprehensive 
model  to  develop  some  scaling  guidance  for  the  design  of  PRS  loads  in  DBCADB. 
.Several  issues,  such  as:  the  low  overall  system  efficiency,  the  possibility  of  POS 
restrike  and  other  MTTL  losses,  and  PRS  peak  voltage,  have  been  addressed. 

The  "gas  bag"  PRS  model  is  used  to  assess  the  peak  kinetic  energy  transferred  to 
the  model  load,  then  a  kinetic  energy  calibration  model,  based  upon  both  the  absolute 
load  mass  and  the  specific  kinetic  energy  (per  ion)  delivered  during  the  rundown 
phase,  is  used  to  estimate  the  expected  yield.  If  Kpeak  is  the  largest  kinetic  energy 
transferred  to  the  load,  and  M"  is  the  total  number  of  ions  in  the  load,  then  the  kinetic 
energy  per  ion  can  be  normalized  to  the  energy  required  to  achieve  K  shell  ionization 
and  emission.  The  resulting  ratio  is  a  measure  of  implosion  quality,  viz.  if 


_  Kpeak  v 

“  jVl.46.E3.5i  >  ^  f 

is  larger  than  1,  then  ample  energy  has  been  supplied  to  excite  the  element’s  K  shell. 
The  theory  goes  further  to  relate  the  expected  yield  to  the  parameter  rj  and  the  total 
mass  relative  to  a  mass  break  point  which  defines  the  boundary  between  efficient 
(Yk  oc  AT)  and  inefficient  (Yk  oc  M^)  yield  scaling,  viz. 
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where  0(77)  and  6(77)  are  calibrated  through  fits  to  the  1-D  MHD  solution  data.  Further 
modifications  to  the  yield  scaling  are  folded  in  when  77  »  1  and  when  the  stagnation 
radius  is  made  larger  than  that  seen  in  the  calibrating  MHD  studies,  but  the  result 
can  be  expressed  as  a  simple  fimction  YKiji,  fi)  which  is  proportional  to  AT  for  /x  >  1 
and  for  <  1. 

•  Power  Flow  in  DECADE  Point  Designs 

The  detailed  DECADE  point  design  used  here  has  a  full  spedfication^^^  elsewhere. 
An  element  by  element  summary  description  is  shown  here  for  the  full  (parallel)  16 
modxile  line  in  Fig.  1.,  and  a  single  DECADE  module  (SDM)  is  obtained  from  this 
circuit  by  multipl3dng  all  the  impedances  by  the  number  of  modules,  including  the 
impedances  of  the  convolute  since  the  SDM  is  presumed  to  drive  an  short  drcmted 
convolute  and  PRS.  In  the  SDM  version  of  the  model  the  current  entering  the  POS 
element  is  not  divided  by  the  number  of  modules. 

In  a  typical  calculation  the  transfer  capadtor  (TC)  is  charged  to  1800  kV,  a  level 
somewhat  higher  than  the  expected  Mara  charge.  The  overcharge  serves  to  make  up 
for  energy  that  would  continue  to  flow  into  the  TC  as  it  discharges  into  the  POS  and 
the  inductive  storage  region  formed  by  the  PFL  MITL,  The  water  output  line  is  long 
enough  that  only  one  pulse  from  the  TC  can  ever  interact  with  the  model  PRS  load, 
the  transit  time  from  TC  to  load,  back  to  the  TC,  and  then  back  to  the  stack  voltage 
monitor  exceeds  300  ns.  So  long  as  model  PRS  nmdown  times  are  less  than  200  ns, 
the  model  line  will  be  an  accurate  representation  of  the  full  piilser  shot,  to  indude 
the  voltage  reversal  and  recoil  of  eneigy  back  to  the  TC  that  occurs  at  about  400  ns 
after  the  onset  of  the  TC  discharge.  The  stack  voltage  monitor  signal  point,  shown  on 
the  figure,  is  a  routine  model  output,  as  is  the  current  shunted  by  the  POS,  and  the 
load  current  at  the  feed  MITI/convolute  boimdary.  Feed  MITL  losses  are  modeled 
with  additional  shunts  and  attenuators  in  the  area  indicated.  The  convolute  and  load 
model  has  two  design  limits  ->  dubbed  "low  risk"  and  "sporty"!^!.  For  transmission 
line  calculations  the  sum  of  the  delay  times  shown  (1.76  ns)  and  the  total  inductance 
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of  the  front  end  options  shown  (102.4  — ♦  71.56  nH)  can  be  used  to  set  an  impedance 
Z  =  X/r  for  a  single  combined  element  lying  between  the  feed  MITL  and  the  PRS 
load  proper.  Widening  this  inductance  range  somewhat  (105  -»  70  nH)  defines  the 
"low  risk"  to  "sporty"  spectrum  of  inductances  used  in  the  energy  coupling  studies. 

Turning  to  the  POS  element,  it  conforms  to  the  PI  design  in  physical  dimension 
and  appears  to  the  primary  transmission  line  as: 

i  an  input  impedance  (which  varies  with  the  time  development  of  the  conduction 
phase  A-K  gap  upstream  of  the  Hall  snowplow), 

ii  a  lumped  shtmt,  and 

iii  an  output  impedance  (which  corresponds  in  a  full  line  model  to  a  16  element 
parallel  combination  of  lines  with  the  listed  anode  and  cathode  radii). 

In  calibrating  this  POS  model  only  one  free  parameter,  the  initial  ion  density,  is 
available.  A  value  of  6.25*  10^^  ions/cm^  when  used  in  the  Hall  snowplow  conduction 
phase  model  established  both  the  proper  conduction  time  (»  290  ns)  and  the  proper 
peak  conduction  current  (»  2.05  MA)  for  the  fixed  length  and  radii  of  the  PI  design. 
This  operating  point  is  quite  insensitive  to  the  subscale  gridding  of  the  POS  model,  a 
21  — » 41  cell  change  in  spatial  resolution  changed  the  peak  conduction  current  by  only 
20  kA.  The  addition  of  the  series  resistance  in  the  MITL  changed  these  current  values 
as  well,  bringing  the  peak  conduction  current  down  to  1.94  MA.  All  these  changes 
are  well  within  the  expected  range  for  enhanced  resolution  of  the  same  physics, 
and,  hence,  we  regard  the  present  operation  of  the  Hall  Snowplow  conduction  phase 
modell^l  as  adequate  for  the  SDM  point  design.  Again,  the  ambient  density  choice  is 
the  only  free  parameter  in  the  simplified  Hall  snowplow  conduction  phase  model,  all 
other  relevant  lengths  and  times  have  been  carefully  set  in  accordance  with  the  PI 
design  specifications. 

•  The  PRS  in  a  DM1  Environment 

The  next  task  in  the  POS/PRS  modeling  effort  was  the  establishment  of  the 
hpaolino  performance  of  a  single  DECADE  module  (SDM)  as  conceived  by  PI.  As 
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discussed  above,  the  initial  prediction  of  POS  conduction  time  and  peak  conduction 
phase  current  was  quite  successful  in  that  a  choice  of  6.25  •  10^^  ions/cc  produced  in 
the  simulation  both  the  correct  conduction  time  (  290  ns)  and  the  correct  conduction 
current  ( 2.0MA). 

Sin(»  the  PRS  is  electrically  isolated  from  the  circuit  during  the  conduction  phase, 
any  study  of  its  properties  should  proceed  from  a  fixed  initial  condition.  In  modeling 
the  DECADE  line,  therefore,  the  POS  conduction  phase  was  computed  up  to  a  time 
just  prior  to  the  opening  event  and  saved  as  a  restart  file.  Having  established  this 
fixed  conduction  phase  database,  the  rundown  and  stagnation  of  a  model  PRS  load  in 
the  SDM  was  computed  over  a  domain  of  several  load  masses  on  the  order  of  25 
200  ng/cm.  Early  runs  showed  a  poor  energy  coupling  to  these  loads,  on  the  order  of  5 
%  of  the  energy  available  in  the  lES  region,  but  this  was  due  to  an  incorrect  input  to 
the  models  front  end  inductance.  The  addition  of  more  inductance  to  the  front  end  in 
modeling  the  SDM  resulted,  as  expected,  in  the  coupling  of  additional  energy  to  the 
load,  whether  it  was  a  matched  load  or  a  model  PRS.  For  example,  taking  Pis  general 
range  of  3 — 7  nH  for  the  parallel  hookup,  choose  5  nH  and  thus  obtain  an  80  nH 
single  element  fixed  load  inductance  (prior  to  the  PRS).  The  energy  coupling  rises  to 
about  35  %  of  the  energy  delivered  to  the  lES  in  the  conduction  phase.  Of  course  the 
lumped  inductance  of  the  front  end  is  only  part  of  the  line  specification.  The  actual 
impedance  jump  at  the  transition  into  the  front  end  is  important  in  establishing  the 
best  model  because  the  reflection  coefficient  for  that  interface  plays  a  role  in  the 
complicated  power  coupling  process  that  follows  the  POS  opening.  Since  PI  furnished 
the  more  detailed  treatment  of  the  SDM  front  end  discussed  above,  allowing  both 
the  Imnped  inductance  and  the  transition  impedance  to  be  taken  into  account,  this 
coupled  energy  fraction  has  remained  higher.  In  all  the  studies  below  which  vary 
the  POS  and  front  end  parameters  into  a  fixed  SDM  PRS  load,  the  overall  front  end 
inductance  was  constrained  to  the  range  of  70  105  nH. 

Having  established  a  reasonable  operating  range  for  the  SDM  line  model,  in  terms 
of  POS  and  PRS  feed  parameters,  the  next  logical  step  was  to  examine  the  details  of 
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the  model’s  operation  for  fixed  POS  parameters  and  a  PRS  mass  near  the  optimum 
for  energy  transfer  to  the  load.  Note  that  optimum  load  cavity  (or  front  end)  energy 
transfer  mass  is  usually  not  the  optimum  mass  for  kinetic  energy  transfer,  but  the 
two  masses  are  close.  Examples  of  the  SDM  operational  scenario  are  shown  below. 
In  "SDM  Line  Profiles"  (Fig.  2)  the  line’s  three  propagation  media  (water,  oil,  and 
vacuum)  are  shown  in  the  "bar  legend"  and  the  abscissa  is  in  time  delay  [ns]  along 
the  transmission  line.  The  current  profile  peaks  at  2.05  MA  just  upstream  of  the 
POS,  which  is  caught  in  the  act  of  opening.  About  half  the  current  resident  in  the 
lES  is  beginning  to  move  toward  the  load.  The  voltage  profile  peaks  at  11.8  MV  just 
downstream  of  the  POS,  but  this  voltage  is  a  factor  of  three  or  more  greater  than 
either  the  peak  load  voltage  or  the  peak  stack  voltage  seen  later.  The  line  impedance 
profile  peaks  at  22.4  Q  at  the  oil/vacuum  transition  and  constrains  the  profile  of 
.  energy  density  to  peak  immediately  upstream  of  the  POS,  viz.  within  the  lES. 

In  SDM  line  Waveforms  (Fig.  3)  are  shown  the  time  series  for  energies,  voltages, 
and  currents.  No  energy  refiects  from  the  matched  load  a  resistively  terminated 
inductance  that  matches  the  local  output  impedance  of  4.5%o  at  the  end  of  the  MTTL. 
The  load  absorbs  about  35  %  of  the  eneigy  delivered  to  the  IBS  by  the  transfer 
capacitor.  The  load  energy  history  and  line  input  energy  history  are  shown  in  the 
red  and  green  curves.  The  load  current  peaks  at  about  2/3  MA,  mostly  due  to  the 
somewhat  larger  inductance  of  the  front  end  parameters.  Note  that  the  load  voltage 
peak  is  more  than  twice  the  stack  voltage  peak,  due  to  different  local  impedances. 

The  surface  plots  compare  the  matched  load  and  PRS  with  respect  to  current, 
voltage,  and  POS  restrike.  In  "Line  Current  In  Space  and  Time"  (Fig.  4),  the  spatial 
scale  is  now  physical  path  length  (  8m),  and  the  broad  yellow  line  cuts  the  suiface  at 
the  oil/vacuum  transition.  The  only  significant  feature  distinguishing  the  two  current 
histories  is  the  null  in  line  current  that  occurs  for  the  matched  load  but  not  for  the 
PRS.  Comparing  the  "Line  Voltage  in  Space  and  Time"  (Fig.  5)  histories,  the  PRS 
produces  slightly  higher  peak  voltages  at  opening  all  along  the  MTTL  section  of  the 
line.  The  MITL  section  is  marked  by  a  more  open  grid  on  the  surface  plot,  arising 


from  the  longer  spatial  increment  for  a  fixed  time  delay.  For  either  load  the  POS 
opening  event  sends  a  reflected  voltage  pulse  back  to  the  TC,  then  to  be  reflected 
again  back  toward  the  load.  The  secondary  reflection  might  be  artificial,  due  to  an 
impedance  jump  in  the  model  where  it  couples  to  the  TC,  but  it  has  no  effect  on  the 
PRS  dynamics  calculation  because  enough  delay  is  built  into  the  model  output  line  to 
provide  isolation.  Comparing  next  the  "Current  Fraction  Shunted"  (Fig.  6)  histories, 
the  existence  of  a  current  null  for  the  matched  load  forces  a  POS  restrike  (at  very 
low  current  of  course)  because  the  low  line  current  cannot  magnetically  insulate  the 
gap,  even  when  the  full  POS  gap  is  quite  large,  as  in  the  PI  design.  The  restrike  here 
shows  the  sensitivity  of  this  particular  numerical  diagnostic  to  the  restrike  condition. 
So,  even  though  the  test  cases  chosen  for  the  PRS  model  tend  to  show  no  serious 
restrike,  the  diagnostic  would  pick  it  up  if  the  model  line  voltage  and  current  solution 
predicted  it.  It  would  appear  that  current  nulls  are  more  dangerous  than  higher 
voltages  in  producing  a  restrike  situation. 

The  next  step  after  matched  load  studies  was  the  addition  of  a  PRS  load  in  a 
mass  range  determined  to  provide  load  stagnation  before  500  ns.  Energy  transfer 
to  the  load  continued  to  be  about  what  was  expected  from  earlier  work.  In  the  case 
of  the  PRS  load  the  details  of  the  "gas  bag"  stagnation  dynamics  are  shown  on  Fig. 
7.  Compared  to  the  matched  load,  for  the  PRS  the  load  current  is  a  bit  larger  (  3/4 
MA),  and  the  energy  absorbed  is  down  to  about  10  %  of  the  peak  transferred  energy. 
Stack  and  load  voltages  are  comparable  to  the  matched  load  case.  The  peak  kinetic 
energy  delivered  to  the  load  is  l.SkJ  ,  while  the  inductance  change  is  13.6nH  at 
peak  compression,  see  below.  A  graphic  summary  of  the  SDM  performance  with  a 
nominal  80nH  front  end  and  a  PRS  load  is  shown  below  (Fig.  8).  In  the  line  plots 
the  normalizing  value  for  any  particular  (color  coded)  quantity  is  shown  in  the  legend 
accompanying  the  graph. 

The  overall  electrical  performance  of  the  SDM  calciilation  appears  to  be  in  fine 
agreement  with  the  calculations  done  at  PI  when  compared  on  the  basis  of  line  current 
and  voltage  waveforms.  As  a  continuing  effort  to  evaluate  the  model  we  furnished  PI 
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a  detailed  graphic  summary  of  the  SDM  model  electrical  output  at  the  relevant  test 
points  for  comparisons.  Once  the  fundamental  electrical  performance  of  the  model 
was  established,  the  examination  of  PRS  load  kinetic  energy  transfer  commenced 
with  a  finer  resolution  scan  of  mass  to  determine  the  optimiun  kinetic  energy  mass 
loading  for  the  inductance  and  impedance  parameters  fixed  by  the  SDM  line  design. 

Several  studies  were  carried  out  to  understand  the  interaction  of  a  POS  with 
a  given  PRS  load.  Operating  near  the  optimum  kinetic  energy  transfer,  one  may 
construct  any  of  several  POS  parameter  scans  to  determine  the  sensitivity  of  the  final 
PRS  kinetic  energy  (T/)  to  POS  behavior.  There  are  three  particular  variations  of 
interest: 

(  i)  Since  density  changes  the  conduction  time,  fix  the  POS  sheath  gap  (Dpos)  and 
vary  the  initial  POS  density  (noc)  to  vary  T/  (n<,c). 

(  ii)  Since  the  equilibrium  sheath  gap  size  (Dpo5)  at  the  start  of  opening  partially 
determines  the  gap  impedance  increase,  fix  the  conduction  time  (through  n^c)  and 
vary  the  gap  to  determine  T/  (Dpos). 

(  iii)  Examine  the  sensitivity  of  the  energy  delivery  to  the  late  time  impedance  of 
the  POS  model. 

As  discussed  above  the  restrike  tendency  (item  iii.),  as  measured  by  the  open  phase 
energy  dissipation  in  the  switch  (Epos),  may  be  larger  for  larger  kinetic  energy 
delivery.  So  one  means  of  investigation  is  to  examine  progressively  lower  mass  cases, 
\mtil  the  motional  impedance  on  the  rundown  can  be  made  large  enoxigh  to  cause 
some  POS  dissipation. 

The  second  SDM  study  reveals  a  clear  tendency  for  the  PRS  to  perform  better  for 
larger  POS  gaps  at  the  start  of  the  opening  phase.  The  initial  gap  is  always  pushed  to 
the  largest  allowed  value  in  the  typical  POS/PRS  calcnalation,  so  the  variation  of  the 
initial  gap  size  might  be  expected  to  vary  only  the  initial  energy  transfer  to  the  load. 
It  turns  out  however  that  the  DM1  front  end  parameters  are  such  that  a  large  value 
for  Dpos  tends  to  produce  a  large  amplitude  "sloshing"  of  energy  into  and  back  out  of 
the  PRS  frond  end.  So  much  so  that  the  PRS  actually  begins  a  dear  implosion,  slows 
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down  as  the  current  leaves  it,  then  surges  on  as  the  next  crest  returns  from  reflection 
near  the  POS.  In  a  well  optimized  rundown  there  are  as  many  as  three  of  these 
sloshing  episodes.  Lowering  Dpos  toward  the  0.5  mm  range  inferred  from  Hawk 
opening  phase  behavior,  one  finds  that  the  energy  sloshing  is  reduced  in  amplitude 
and  slightly  lengthened  in  period.  The  peak  energy  delivered  to  the  PRS  front  end, 
and  the  peak  kinetic  energy  of  the  PRS,  and  the  peak  load  current  decay  modestly  in 
the  process.  The  behavior  is  shown  in  Fig.  9.  Notice  that  the  peak  load  current  in 
these  transients  can  exceed  the  POS  conduction  phase  peak  current. 

Completing  the  program  of  study  just  laid  out,  the  first  of  the  SDM  studies  cited 
above  (load  kinetic  energy  with  POS  density)  reveals  a  clear  tendency  for  the  PRS  to 
perform  better  for  early  "firing  angles",  as  determined  by  the  initial  POS  density.  As 
the  conduction  phase  is  made  shorter  and  the  POS  made  to  open  at  smaller  current 
level,  the  kinetic  energy  to  the  load  rises  by  about  20  percent,  as  shown  in  Fig.  10. 
These  studies  only  varied  the  POS  density  and  that  density  cannot  be  taken  to  zero. 
The  strict  direct  drive  case  thus  required  a  slight  code  modification  and  it  was  added 
in  the  full  line  simulation,  allowing  a  true  direct  drive  series  to  be  done  as  reported 
below. 

•The  PRS  in  a  Full  DECADE  Environment 

Combining  16  DM1  modules  in  parallel  to  form  the  full  DECADE  machine  model  is 
a  simple  matter  of  dividing  the  line  impedances  by  Nmodu/o  the  number  of  such  units 
used.  In  addition,  the  front  end  inductance  is  divided  by  Nmodttie  while  keeping  the 
same  time  delay,  in  order  to  transform  to  the  parallel  impedance.  The  current  feeding 
the  POS  model  is  also  divided  by  limoduu  to  preserve  the  conduction  time  and  peak 
current  in  each  POS  module.  On  the  POS  output  side,  the  voltage  across  the  shunts 
is  preserved  by  maintaining  the  same  relative  voltage  decay  factor  that  occurs  for  a 
single  module  POS  model  in  equivalent  electrical  circumstances.  Elsewhere  in  the 
circuit,  the  transfer  capacitor  (TC)  is  transformed  to  maintain  the  same  quarter  (yde 
time,  and  each  line  shunt  or  attenuator  is  also  adjusted  to  the  parallel  equivalent. 
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The  net  effect  of  the  transformation  is  a  single  line  model  carrying  Nmotfuie  times 
more  current  into  the  array  of  POS  units.  When  the  POS  modules  open,  their  parallel 
impedance  determines  the  energy  dissipated  and  transfered  downstream  to  the  PRS 
load  that  now  experiences  the  implosion  force  generated  by  the  combined  currents. 
Similar  to  the  behavior  seen  in  modeling  DM1,  the  POS  opening  event  in  DECADE 
creates  qmte  a  bit  of  ringing  in  the  front  end  between  the  POS  and  the  load,  and 
the  transient  excursions  of  current  about  the  low  frequency  component  can  lead  to 
load  currents  in  the  30  MA  range.  The  high  transient  peak  current  is  smoothly 
degraded  by  the  rundown  of  the  PRS  however,  and  typical  drive  currents  just  prior 
to  stagnation  are  in  the  12  — »  18  MA  range.  There  is  also  some  sensitivity  to  the 
early  resistance  value  inferred  for  the  Pf^  sltig  phase,  a  factor  10  change  in  this 
resistance  can  lead  to  a  factor  of  almost  2  in  the  expected  K  shell  yield  for  loads  close 
to  the  *7  >1  domain.  However,  as  the  kinetic  energy  rises  into  the  domain  sufficient 
for  substantial  K-shell  yield  this  sensitivity  gets  weaker. 

The  useful  comparison  of  two  point  designs  differing  only  in  front  end  inductance 
on  the  load  side  of  the  convolute  can  be  established  by  examining  a  similar  array  of 
load  masses  for  each  design.  The  POS  behaviour  is  fixed  in  both  cases  by  calculating 
from  a  restart  file  taken  at  the  end  on  the  conduction  phase  and  simply  shifting  the 
downstream  load  and  convolute  model  components. 

For  a  lower  inductance  front  end  the  study  reveals  that  it  takes  (roughly)  a  40  % 
drop  in  the  inductance  to  achieve  a  8  •  10  %  increase  in  the  yield,  for  the  higher  mass 
cases.  Conversely,  any  convolute  inductanro  in  the  proper  range  can  therefore  be 
increased  somewhat  if,  for  example,  that  increase  were  designed  around  allowing  a 
longer  implosion  distance  for  the  PRS  -  a  modification  which  could  improve  the  yield 
by  a  much  larger  percentage. 

Moreover,  the  peak  kinetic  energy  transfer  is  found  to  occur  for  a  lower  mass  as 
the  convolute  inductance  is  increased.  A  higher  inductance  implies  a  lower  kinetic 
energy  at  any  fixed  mass,  while  a  lower  inductance  gives  a  better  compression  for  the 
pinch  and  a  slightly  larger  change  in  load  inductance  prior  to  stagnation. 
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•  The  PRS  in  a  Direct  Drive  DECADE  Environment 

Unlike  the  above  study  of  PRS  and  POS  performance  as  the  POS  density  was 
decreased  to  change  the  "firing  angle"  of  the  DECADE  line,  a  true  direct  drive  removes 
the  POS  completely,  viz.  there  is  no  conduction  phase  and  no  attempt  to  sharpen 
the  power  delivery  from  the  transfer  capacitor.  This  is  a  very  attractive  option  from 
the  standpoint  of  operational  simplicity,  because  essentiaUy  all  the  electrical  and 
mechanical  measures  of  PRS  performance  are  improved  at  higher  mass  loadings  over 
even  the  best  case  scenario  based  on  the  removal  of  front  end  inductance.  At  lower 
mass  loadings  the  comparison  swings  back  to  favor  the  POS  configuration  slightly, 
but  in  no  case  does  the  direct  drive  option  exact  a  m^or  penalty  on  kinetic  energy 
delivery. 

The  comparison  can  be  refined  further  by  examining  the  transfer  of  kinetic  energy 
in  the  two  options.  In  Fig.  11  are  shown  contour  plots  of  the  peak  kinetic  energy 
(Kptak)  achieved  in  the  two  cases.  The  POS  option  offers  a  more  sustained  delivery  for 
lower  masses,  due  to  its  quick  opening  time,  on  a  wide  domain  of  initial  radii;  while 
the  direct  drive  option  offers  a  more  regular  variation  in  delivered  energy,  due  to  the 
absence  of  reflections  in  the  front  end  between  the  POS  and  the  convolute.  The  same 
reflections  were  responsible  for  the  surge  behavior  seen  in  the  SDM  studies  discussed 
above.  The  POS  option  also  offers  a  modest  increase  in  delivered  kinetic  energy  on  the 
domain  [mR^  »  200,  R  »  5],  but  this  performance  variation  is  undoubtedly  contingent 
on  the  accuracy  of  the  POS  model  in  terms  of  opening  time  and  voltage  standoff. 

When  the  yield  scaling  YKiv^  is  examined  for  these  two  kinetic  energy  delivery 
options,  the  resulting  yields  are  similar  for  low  Z  materials,  c.f.  Fig.l2(a,b),  and 
marginally  favorable  to  the  POS  driver  for  higher  Z  materials,  c.f.  Fig.l2(c,d).  Given 
the  similarities  in  yield  performance  shown  here  and  the  level  of  imcerte^nt  in  POS 
operation,  a  DECADE  line  built  to  drive  A1  or  Ar  would  almost  certainly  be  better 
configured  as  a  direct  drive  system.  If  the  primary  load  domain  is  moved  up  to  Cu 
or  Kr,  then  there  is  clearly  something  to  be  gained  in  POS  mode.  The  behavior  of 
the  POS  driver  with  respect  to  restrike,  redosure,  voltage  holdoff  and  overall  energy 
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dissipation  then  becomes  the  pivotal  question.  A  POS  placed  too  close  to  the  load 
convolute  could  lose  all  the  advantage  shown  in  these  figures  by  injecting  plasma  into 
the  feeds  and  either  slowing  or  diminishing  the  energy  flow  to  the  PRS. 

•  Project  Guidance 

The  present  studies  have  demonstrated  that  DECADE  is  reasonably  sure  of 
reaching  the  early  PRS  radiation  goals  assigned  to  it,  e.g.  200  kJ  of  Cu  K  shell. 
The  POS  is  of  at  least  marginal  utility  in  achieving  this  end,  provided  that  it  does 
not  expend  much  plasma  into  the  downstream  MITL  and  convolute.  Should  the 
"sneezing"  POS  present  such  a  problem,  a  simple  fix  would  be  to  operate  in  direct 
drive  mode,  with  only  a  modest  drop  in  the  yield  for  loads  like  Cu  and  a  more  serious 
reduction  for  Kr.  A  more  complete  fix  would  involve  a  focus  on  two  follow  up  study 
.  areas  with  a  more  detailed  line  model. 

The  first  area  would  be  the  investigation  of  kinetic  energy  transfer  and  convolute 
impedance  changes  corresponding  to  modifications  in  the  final  PRS  feeds.  The  focus 
would  be  to  demonstrate  that  a  wider  stalk,  and  the  larger  initial  load  radius  that 
would  allow,  could  be  tolerated  with  regard  to  the  added  inductance  very  close  to 
the  load.  The  experience  to  date  suggests  that  such  a  wider  stalk  would  be  a  small 
perturbation  to  the  overall  power  flow  and  a  great  help  in  getting  better  kinetic 
energy  transfer  to  the  pinch.  In  addition  the  examination  of  axial  magnetic  fields  to 
stabilize  the  long  rundown  options  would  clearly  be  advisable. 

The  second  area  involves  the  improvement  of  the  POS  model  to  accept  a  new  gap 
dynamics  equation  (developed  this  year,  but  not  tested  extensively).  The  POS  may 
be  placed  closer  to  the  convolute  in  order  to  shorten  the  risetime  into  the  load,  and 
the  new  operating  domain  shoiild  be  investigated  with  regard  to  the  possibility  that 
the  POS  will  physically  redose  due  to  the  stresses  of  the  larger  applied  voltages.  In 
addition  a  more  complete  MITL  loss  model  could  be  added  to  assess  the  impact  of 
downstream  POS  material. 

The  knowledge  developed  through  these  additional  studies  would  also  help  in 


59 


much  needed  studies  supporting  the  Jupiter  concept. 
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Fig.2  SDM  Line  Profiles  in  V,l  at  290[n8],  just  as  the  POS  opens. 
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Fig.  3  SDM  Line  Waveforms  on  [0  -  0.500  ^s]  with  a  matched  load. 
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Fig.  9  POS  gap  size  and  KE  performance 

*  Large  gaps  at  opening  set  up  surges  to  the  load 
and  result  in  a  better  KR  transfer. 
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Fig.  12  Yields  from  Direct  Drive  and  POS  modes  in  Decaoe 
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V.  ANALYSIS  OF  MHD  DATA 

Predictions  of  DECADE  performance  and  the  design  of  better  loads  to  improve  performance 
depend  on  the  ability  to  successfully  match  theoretical  models  with  experimental  data.  Besides 
predicting  gross  yields  and  implosion  times  and  their  relation  to  shot  parameters,  theoretical  models 
should  accurately  follow  the  load  dynamics  during  the  pinch.  More  detailed  comparisons  between 
machine  performance  and  model  predictions  would  make  it  easier  to  identify  where  theoretical 
models  should  be  extended  or  where  load  characteristics  should  be  modified.  This  will  increase 
reliability  and  lead  to  more  acurate  theoretical  models. 

Even  in  the  absence  of  detailed  experimental  measurements,  the  complexity  of  numerical 
modeling  makes  it  essential  to  develop  the  highest  quality  analytical  tools.  Analysis  of  MHD  data 
is  very  similar  to  that  of  experimental  data:  Both  numerical  models  and  experimental  instruments 
produce  i^w  data  that  has  little  direct  or  obvious  physical  significance.  The  meaning  in  terms 
of  physical  quantities,  and  underlying  processes  must  be  deduced  by  analyzing  this  raw  data  and 
making  inferences.  The  challenge  in  analyzing  experimental  data  is  to  deduce  plasma  conditions 
from  measurements  like  spectra  and  fields,  while  the  challenge  in  analyzing  MHD  data  is  to  infe'- 
experimental  measurements  from  numerically  determined  plasma  conditions. 

A  recent  series  of  experiments  by  G.  Davara,  M.E.  Foord,  Y.  Maron  and  collaborators  at  the 
Weizmann  Institute  of  Science^  has  utilized  both  radially  and  axially  resolved  measurements  of 
spectra  and  the  process  of  ionization  to  obtain  information  about  heating,  compression  and  transport 
within  the  pinch.  Such  spatially  resolved  measurements  (especially  the  radial  profiles)  have  rarely 
if  ever  been  obtained  before  for  a  z-pinch  implosion,  and  provide  useful  information  about  pinch 
dynamics.  The  quantities  measured  are  easily  computed  using  numerical  models,  and  comparison 
of  such  measurements  with  model  predictions  could  be  an  important  check  on  the  models. 

Our  PRS  model  includes  full  1-D  MHD  mass  transport  as  well  as  detailed  ionization  and 
radiation  dynamics,  and  is  thus  highly  nonlinear.  Small  deviations  of  pinch  parameters  from 
model  predictions  could  result  in  large  discrepancies  between  predicted  and  actual  implosion  times 
or  radiative  yields,  both  because  the  dynamics  are  nonlinear  and  because  new  processes  come 
into  play.  For  example,  viscosity  is  dependent  on  plasma  conditions,  the  electron  distribution 
and  the  presence  or  absence  of  turbulence;  different  viscosities  will  result  in  different  implosion 
times,  different  temperatures  and  degrees  of  compression,  and  so,  in  a  ijonlinear  way,  in  different 
yields  and  spectra  of  radiation.  The  radiative  component  of  this  is  particularly  important,  because 
the  pinch  radiates  strongly,  particularly  at  stagnation,  and  because  the  radiative  output  depends 
sensitively  on  the  degree  of  ionization  attained  in  the  plasma,  which  in  turn  depends  on  the  peak 
plasma  conditions.  Radiation  is  also  affected  by  the  electron  distribution  function  in  the  interior, 
and  this  has  been  shown  to  be  critically  dependent  on  the  conductivity  and  degree  of  magnetization, 
which  in  turn  depend  on  the  mobility  of  the  electrons.  Electron  runaways  could  also  result  from 
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rapid  radial  current  penetration,  decoupling  the  load  from  the  generator  and  inducing  hydrodynamic 
instabilities  farther  along  the  axis.  We  have  previously  shown*  that  runaway  electron  production  is 
sensitively  dependent  on  the  electric/magnetic  field  ratio. 

In  this  report,  we  will  first  briefly  describe  the  results  of  the  Weizmann  Institute  experiment  and 
then  show  the  analogous  results  obtained  from  a  detailed  ID  MHD-based  radiative-hydrodynamic 
modeP  of  an  implosion  with  parameters  appropriate  for  DECADE.  The  presentation  is  not  meant 
to  be  a  direct  comparison  between  the  two  sets  of  data.  This  would  not  be  appropriate,  because  of 
the  great  difference  in  implosion  parameters  in  the  two  cases.  Rather,  it  is  meant  to  demonstrate 
the  capability  of  our  numerical  model  to  computationally  obtain  the  same  kind  of  data  as  the 
experiment.  These  data  could  serve  as  a  point  of  comparison  between  model  and  experiment,  if 
similar  measurements  could  be  made  on  high-power  machines  as  were  made  on  the  Weizmann 
apparatus.  The  data  are  also  useful  in  themselves,  because  they  give  insight  into  the  complex 
ionization  dynamics  described  by  the  model.  Thus,  the  analysis  is  part  of  a  continuing  program  to 
assure  self-consistency  and  to  obtain  maximum  physical  insight  from  the  model  results.  This  kind 
of  analysis  will  be  extended  in  the  future  to  provide  a  more  intuitive  understanding  of  implosion 
dynamics  in  DECADE,  Saturn  and  other  machines. 

Experimental  results 

The  Weizmann  Institute  experiments  utilize  a  low  power,  carbcn  dioxide  gas  puff  z-pinch. 
The  load  is  injected  with  initial  radius  R  fa  2  cm  and  annular  extr  .it  Sr  fa  5  mm.  The  peak  load 
current  Imar  ^  250  kA,  with  fa  30  kV  potential  difference,  and  the  current  rise  time  is  about 
1.1  fisec.  The  implosion  timescale  is  about  600  ns. 

The  experimenters  used  a  1.3  m  spectrograph  (resolution  0.04  A)  coupled  to  a  high-resolution 
streak  camera  (1(X)  ps  temporal  resolution),  with  sensitivity  between  2000  A  and  8000  A.  The 
spatial  resolution  was  0.5  mm  horizontally  and  2.0  nun  vertically.  Relevant  to  the  present  report, 
they  measured  the  radial  position  and  ocity  of  oxygen  charge  states  in  the  load  (O  n  through 
O  VI).  The  time-evolution  of  the  charge-state  position  showed  the  speed  of  ionization,  while  the 
Doppler  broadening  in  the  spectrum  showed  the  velocity  of  the  ions. 

The  experimenters  found  that  each  charge  state  behaved  individually.  The  ionization  front, 
defined  by  the  2^}pearance  of  O  D,  advanced  much  more  riq)idly  than  any  Doppler  shift-derived 
species  velocity,  and  reached  a  maximum  of  10  cm//isec  just  before  stagnation.  Initially  neutral, 
the  gas  progressively  ionized  at  the  outer  edge,  which  was  the  hottest  point  (18  eV  at  80  ns  before 
stagnation).  Species  velocity  was  much  faster  for  ions  at  higher  degrees  of  ionization;  the  O  11 
fnaximum  velocity  was  about  4  cm//tsec,  while  O  VI  reached  a  maximum  of  10  cm//isec  just  before 
stagnation. 

The  features  that  will  coiKxm  us  here  are  an  ionization  wave  propagating  faster  than  the  local 
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fluid  velcKity,  and  the  faster  kinetic  velocity  of  higher  charged  states  (probably  related  to  the 
localization  of  ionization  at  the  outer  edge).  Ionization  was  surmised  by  the  Weizmann  Institute 
group  to  be  driven  by  electron  heat  conduction,  and  the  ion  fluid  velocity  derived  primarily  from 
ion-ion  collisions  rather  than  a  pressure  gradient.  Our  model  results  for  a  high-power  aluminum 
implosion  behave  similarly  to  the  Weizmann  Institute  COj  implosion,  although  the  story  is  more 
complicated.  The  ionization  front  velocity  exceeds  the  fluid  velocity,  and  the  higher  charge  states 
have  greater  velocity  than  lower  charge  states,  but  in  the  model  results  not  all  ionization  occurs  at 
the  outer  edge. 

Model  results 

Here,  we  will  show  calculations  of  radially-resolved  ionization  state  populations  and  plasma 
conditions  from  a  numerical  (ID  MHD)  model  of  a  high-power  implosion  of  aluminum.  These 
numerical  results  are  the  same  quantities  reported  from  the  experiments,  “measured”  to  arbitrary 
precision  with  the  theoretical  model. 

The  1 D  MHD  model  used^  has  been  used  extensively  to  model  implosions  from  a  variety  of  PRS 
devices.  Besides  MHD  particle  dynamics,  it  includes  full  radiation  and  ionization  hydrodynamics; 
for  these  calculations,  55  atomic  states  were  included  in  the  model,  including  all  14  possible 
ionized  states  of  aluminum  and  many  of  their  possible  excited  states.  A  lumped  circuit  model  was 
used,  in  this  case  with  Saturn-based  circuit  parameters.  The  model  utilizes  enhanced  transport, 
so  that  transport  coefficients  have  been  phenomenologically  adjusted  to  match  observed  implosion 
parameters:  electrical  resistivity  is  set  to  20  times,  thermal  resistivity  is  set  to  30  times,  and  viscosity 
is  set  to  80  times  the  Spitzer  value. 

The  load  was  an  aluminum  plasma  with  mass  mass  loading  of  330  /ig/cm,  length  2  cm,  and 
radius  between  1.18  and  1.23  cm  (Sr  =  0.5  mm).  Initially,  the  temperature  is  about  10  eV,  the  ion 
density  is  1 .8  x  10^®  cm"^,  and  the  plasma  consists  only  of  A1 V  and  A1 VI. 

A  plot  of  the  load  current  and  voltage  during  the  modelled  implosion  is  shown  in  Fig.  (1).  The 
time  f  =  0  is  defined  to  be  when  the  voltage  is  initially  applied.  The  current  reaches  a  maximum  of 
8.8  MA  at  69  ns,  and  was  arbitrarily  cut  off  at  the  point  of  maximum  compression.  The  maximum 
voltage  across  the  load  is  4.6  MV,  at  50  ns.  The  average  ion  density  and  temperature  during  the 
pinch  are  shown  in  Fig.  (2).  These  reach  a  peak  nearly  together  at  about  81  ns,  which  is  also 
the  time  of  peak  compression  (stagnation).  At  this  time,  the  pinch  outer  radius  is  0.28  cm.  At 
stagnation,  the  local  electron  density  reaches  1. 1  x  the  electron  temperature  at  the  inner  edge 
is  near  10  keV,  and  the  plasma  consists  mostly  of  completely  stripped  aluminum. 

The  progress  of  the  implosion  can  be  seen  graphically  in  Fig.  (3).  Here  an  RGB  palette  is  used, 
where  black  is  zero,  red  is  the  minimum  finite  value  and  blue  is  the  maximum  value.  The  plots  in 
this  figure  are  analogous  to  streak  photogriq>hs,  with  the  ordinate  corresponding  to  time  and  the 


78 


abscissa  to  radial  position  (the  center  line  is  at  the  origin).  In  all  plots,  the  overall  compression  of 
the  plasma  to  a  minimum  at  80  ns  is  evident.  The  temperature  plots  show  the  ohmic  heating  of  the 
plasma  initially  at  the  outer  edge  and  later  compressional  heating  at  the  center.  The  density  plot 
shows  that  the  pinch  initially  expands  somewhat  and  cools  and  then,  during  run-in,  experiences 
strong  compression  at  the  outer  edge.  The  degree  of  ionization  plot  shows  the  extremely  rapid 
propagation  of  the  ionization  wave  soon  after  compression  begins.  The  electric  and  magnetic  field 
plots  show  the  current  penetration  during  run-in  and  stagnation. 

For  computational  purposes,  the  load  is  divided  into  30  concentric  annular  zones,  each  of 
which  moves  with  the  local  fluid  velocity,  each  constrained  by  the  motion  of  the  others.  Because 
of  this,  the  local  fluid  velocity  is  always  known  and  can  be  used  as  a  benchmark.  If  the  ionization 
wave  propagates  towards  inner  zones  (as  it  does),  its  velocity  exceeds  the  local  fluid  velocity.  The 
outer  radius  of  the  outer  zone  is  taken  to  be  the  pinch  outer  radius,  and  is  close  to  what  a  pinhole 
photograph  would  measure  as  the  pinch  radius,  although  the  luminosity  of  each  zone  has  not  been 
included.  The  motion  of  the  zone  boundaries  during  the  implosion  is  shown  in  Fig.  (4).  Where  the 
boundaries  come  closer,  the  local  density  is  larger. 

The  color  plots  in  Fig.  (5)  show  the  process  of  ionization  in  the  load.  Here  again,  the  ordinate 
is  the  time  axis,  with  a  full  extent  of  100  ns  lasting  until  just  slightly  after  stagnation.  The  abscissa 
corresponds  to  the  radial  zone  number,  and  is  divided  into  30  strips,  one  for  each  zone.  The  pinch 
center  line  is  at  the  origin,  and  higher  zone  numbers  are  closer  to  the  outer  edge.  Using  zone 
instead  of  radial  position  for  the  abscissa  has  two  advantages.  First,  it  expands  the  display,  making 
variation  clearer,  particularly  near  stagnation  when  the  plasma  is  compressed.  Second,  it  shows 
explicitly  that  the  ionization  front 's  moving  in  excess  of  the  fluid  velocity  (that  is,  the  cell  velocity), 
and  by  how  much.  An  RGB  pclette  is  used  here,  and  the  appearance  (and  disappearance)  of  each 
charge  state  is  marked  by  the  appearance  of  color  against  the  black  background.  Colors  represent 
the  magnitude  of  the  population  fraction  of  the  given  charge  state,  and  all  population  fractions  for 
the  different  charge  states  sum  to  one  at  any  given  time  and  any  given  radial  position.  Blue  means 
a  low  population  fraction,  and  red  means  the  plasma  at  that  position  and  time  is  almost  entirely 
composed  of  the  given  charge  state.  Plots  are  aligned  horizontally  so  that  the  plasma  composition 
at  a  given  time  can  be  easily  identified,  at  least  for  ionizated  states  displayed  in  the  same  column. 

It  can  be  seen  in  Fig.  (5)  that  the  plasma  is  initially  composed  of  A1 V  and  A1 VI.  After  about 
30  ns,  the  ionization  wave  commences,  so  that  A1  Vn  appears  at  30  ns,  A1  VIE.  at  35  ns,  A1  DC  at 
40  ns,  A1  X  at  41  ns,  and  A1 XI  at  42  ns,  all  initially  at  the  outer  edge.  Before  about  40  ns,  there 
is  very  little  compression,  and  in  fact  the  plasma  slightly  expands  until  this  time.  After  40  ns,  the 
current  ramps,  and  the  implosion  begins  (this  coincides  with  the  peak  in  the  load  voltage).  The 
increase  in  temperature  and  pressure  during  run-in  rapidly  (within  5  ns)  lead  to  the  development 
of  the  A1 XI  and  A1  Xn  species  at  the  outer  edge,  and  these  species  dominate  until  stagnation  sets 


79 


in,  near  60  ns.  The  Weizmann  Institute  observations  were  made  during  this  run-in  period  in  their 
experiment. 

As  stagnation  occurs,  a  secondary  ionization  wave  propagates  from  the  outer  edge  to  the  center, 
and  A1 XII  further  ionize.*?  to  A1  Xni;  after  about  one  ns,  another  ionization  front  forms,  moving  from 
the  center  out  to  the  edge.  This  last  process  of  ionization  occurs  during  the  maximum  compression, 
as  the  implosion  kinetic  energy  thermalizes,  and  it  produces  fully  stripped  ions.  Although  there  is 
some  cooling  on  the  outside  after  stagnation,  the  plasma  is  mostly  fully  stripped  for  the  remainder 
of  the  run  (to  1(X)  ns). 

Quantitative  details  of  the  ionization  dynamics  can  be  seen  in  the  population  profiles  in  Fig.  (6). 
As  seen  in  the  figure,  population  fractions  at  21  ns  are  essentially  as  in  the  initial  state.  At  30  ns,  the 
outer  edge  has  heated  sufficiently  to  produce  significant  populations  of  A1 VI  and  A1 VII.  At  40  ns, 
outer-edge  ionization  continues,  with  charge  states  to  A1 X  and  A1 XI  forming  there.  Also  evident 
at  this  time,  however,  is  interior  pre-ionization,  so  that  near  1  cm  there  is  more  A1  VI  than  A1  V. 
This  is  a  consequence  of  nonuniform  compression  during  run-in,  which  causes  a  slight  pre-heat  of 
the  pinch  interior  and  may  decrease  the  efficiency  of  thermalization  at  stagnation.  The  main  effect, 
however,  is  the  outer-edge  ionization,  which  continues  to  be  evident  at  62  ns  with  the  formation 
there  of  A1  Xn. 

The  average  radial  velocity  of  the  ions  in  each  charge  state  as  a  function  of  time  during  the 
run-in  is  shown  in  Fig.  (7).  This  is  what  would  be  determined  by  measuring  Doppler  broadening 
of  ionized-state  spectra,  as  was  done  in  the  experiment.  As  was  experimentally  found,  the  higher 
charge  states  have  a  higher  velocity.  This  is  not  unexpected,  since  the  each  successive  charge  state 
initially  appears  near  the  outer  edge,  and  the  outer  edge  has  the  highest  radial  velocity  during  run-in 
(this  can  be  inferred  from  the  zone  boundary  trajectories  in  Fig.  (4)).  The  difference  in  velocity 
between  different  charge  states  was  proportionately  much  greater  in  the  experiment  than  in  the 
model  results,  however.  The  absolute  run-in  velocities  in  the  high-power  model  implosion  are  5  to 
10  times  greater  than  the  experimental  values  (up  to  60  cm/fisec). 

The  location  of  the  maximum  population  of  each  charge  state  is  shown  in  Fig.  (8)  as  a  function 
of  time  during  run-in.  This  is  where  the  charge  state  is  most  likely  to  be  located  in  an  experimental 
measurement  (As  Fig.  (5)  shows,  the  charge  states  are  pretty  well  localized  during  run-in).  The 
trajectory  of  the  charge  state  maxima  traces  out  the  progression  of  the  ionization  front  for  each 
stage  of  ionization.  The  closer  in  space  and  time  the  maxima  for  different  charge  states  are,  the 
faster  the  ionization  process.  The  model  results  agree  with  the  experimental  observations,  that  the 
ionization  front  velocity  exceeds  the  fluid  velocity,  and  that  higher  ionization  stages  appear  first  at 
the  outer  edge  of  the  pinch.  In  the  experiment,  however,  the  radial  separation  between  different 
charge  states  was  much  greater;  the  model  shows  a  much  faster  ionization,  more  like  a  shock. 
Besides  the  greater  power  and  faster  pinch  in  the  model  implosion,  a  possible  explanation  for  this 
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is  the  greater  number  of  ionized  states  available  in  aluminum  compared  to  oxygen. 

Conclusion 

We  have  here  demonstrated  the  application  of  several  new  perspectives  for  analyzing  the  results 
of  numerical  models  of  PRS  implosions.  This  analysis  was  prompted  by  the  recent  demonstration 
at  the  Weizmann  Institute  of  new  diagnostic  measurements  in  z-pinch  implosions,  and  the  results 
we  have  obtained  are  in  the  same  pattern  as  those  experimental  results. 

From  the  MHD  data  analysis,  we  find  that,  as  in  the  experiment,  there  is  an  ionization  wave  that 
propagates  faster  than  the  fluid  velocity,  much  faster  during  run-in,  but  probably  not  so  fast  as  to  be 
a  shock  front.  We  also  find  that  the  higher  charge  state  ions  have  a  greater  velocity.  The  difference 
in  velocity  between  different  types  of  ions,  however,  is  much  smaller  in  the  MHD  results  than  in  the 
experiment.  In  addition,  the  MHD  analysis  showed  nonuniformities  within  the  pinch  that  caused 
•  some  internal  preheat  and  pre-ionization,  probably  due  to  compression.  These  nonuniformities 
could  have  been  present  in  the  experiment  as  well,  but  not  observed,  because  of  the  small  number 
of  data  points  available  experimentally. 

If  similar  radially  and  axially  resolved  measurements  were  obtained  from  high-power  PRS 
loads  such  as  DECADE  and  Saturn,  the  procedures  we  have  followed  here  would  allow  a  closer 
examination  of  the  correspondence  between  the  actual  observations  and  the  theoretical  model  of 
the  experiment.  This  would  improve  the  theoretical  models  and  provide  better  understanding  of 
the  experiment.  Even  without  direct  comparison  to  experiment,  however,  this  analysis  is  useful. 

MHD  data  analysis  such  as  that  demonstrated  here  provides  direct  physical  insight  into  the 
actual  working  of  the  model,  and  shows  among  other  things  where  and  for  what  parameters  internal 
heating  is  ocurring,  what  the  important  heating  and  conduction  processes  are,  and  where  the 
important  ionized  states  first  t^)pear.  Data  from  the  model  also  indicate  whether  conditions  are 
suitable  for  other  processes  to  be  important,  such  as  runaways  and  turbulence.  This  process  of 
identifying  the  important  processes  and  conditions  is  just  as  important  in  a  complicated  numerical 
model  such  as  is  needed  for  advanced  PRS  design  as  it  is  in  an  experimental  procedure.  The 
implications  of  this  analysis  of  MHD  data  have  only  been  hinted  at  in  this  presentation,  because 
the  analysis  is  still  in  I*s  initial  stages,  but  this  process  will  be  further  developed  in  the  future. 
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Figure  captions 

1 .  Current  through  the  load  (solid  line)  and  voltage  across  the  load  (dotted  line)  as  a  function  of 
time.  The  time  <  =  0  is  defined  to  be  when  the  voltage  is  first  applied.  The  current  was  arbitrarily 
cut  off  when  the  pinch  reached  maximum  compression. 

2.  Average  ion  density  (solid  line)  and  ion  temperature  (dotted  line)  as  a  function  of  time.  The 
average  is  given  by  summing  the  quantity  for  all  particles  and  then  dividing  by  the  total  number  of 
particles. 

3.  Color  display  of  various  state  variables  during  the  implosion,  as  a  function  of  time  (ordinate) 
and  radial  position  (abscissa).  The  center  line  is  at  the  origin,  and  radial  profiles  at  a  given  time  are 
along  a  vertical  line  intersecting  the  appropriate  time  value  on  the  ordinate.  Black  indicates  zero, 
and  a  “rainbow”  color  palette  is  used,  with  blue,  yellow  and  red  respectively  indicating  minimum, 

<  intermediate  and  maximum  values. 

4.  Position  as  a  function  of  time  of  the  outer  boundaries  of  each  of  the  30  concentric  annular 
zones  in  the  computational  grid.  The  outer  and  inner  zone  boundaries  are  solid  lines,  and  the  other 
boundaries  are  dotted  lines.  Greater  compression  is  indicated  when  the  boundaries  move  closer. 

5.  Color  display  of  the  evolution  of  the  population  fraction  profiles  of  each  ionized  state  as  a  function 
of  time.  The  ordinate  is  the  time  (in  ns),  and  the  abscissa  is  the  zone  number,  corresponding  to  radial 
position  (see  Fig.  (4)),  with  the  center  line  at  the  origin.  Population  fractions  range  from  zero  to 
one,  and  the  sum  of  the  fractions  of  all  chaige  states  is  one  at  every  position  and  time  where  there  is 
plasma.  Black  indicates  the  absence  of  the  chaige  state,  and  a  color  if  present  indicates  the  relative 
abundance  of  the  charge  state,  with  abundance  increasing  from  blue  to  yellow  to  red.  A  profile 
of  the  plasma  composition  at  any  given  time  can  be  seen  along  a  vertical  line  that  intersects  the 
ordinate  at  that  time  value.  Displays  are  horizontally  aligned  to  permit  comparison  of  simultaneous 
populations  (at  least  if  the  charge  states  are  in  the  same  column).  Because  the  abscissa  is  zone 
number  and  not  actual  radius,  these  profiles  show  no  compression  during  the  implosion,  although 
from  other  figures  it  is  clear  that  the  zones  do  get  significantly  closer  during  run-in. 

6.  Population  fractions  for  important  charge  states  as  a  function  of  radius,  at  times  of  21,  30, 41 
and  62  ns.  Charge  state  fractions  not  shown  are  not  significantly  present.  The  population  fraction 
is  plotted  as  a  dashed  line  with  the  charge  state  identification  printed  near  its  highest  point.  For 
example,  at  21  ns,  the  plasma  is  mostly  A1  V  (peak  fraction  0.9  near  the  edge),  with  some  A1  VI 
and  A1 IV.  Population  fractions  at  each  point  sum  to  one. 
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7.  Radial  velocities  (in  cm//isec)  of  ions  in  different  charge  states,  as  a  function  of  time  during 
run-in.  Here,  symbols  for  highest  charge  states  are  joined  by  lines.  These  are:  “+”  for  A1 XIV, 
for  A1  Xm,  for  A1 XH,  diamond  for  A1 XI,  triangle  for  A1 X,  square  for  A1 IX,  “X”  for  A1 VIII, 
and  for  A1  Vn.  In  addition,  the  “+”  symbol  not  joined  by  lines  represents  A1  VI. 

8.  Radial  position  (in  cm)  of  the  peak  population  of  several  chai:ge  states  as  a  function  of  time, 
during  run-in.  The  solid  line  represents  A1 XIV,  which  is  not  significantly  present  during  this  time. 
Other  charge  states  are  represented  by  the  following  symols:  A1  XIII  (“+”),  A1  XII  (“*”),  A1  XI 
(“.”).  A1 X  (diamond),  A1  DC  (triangle),  A1  VIE  (square),  and  A1  Vn  (“X”).  Higher  charge  states 
appear  first  on  the  outer  edge  of  the  pinch,  and  generally  remain  farther  out  from  the  center  line 
during  run-in.  Some  mixing  of  peak  positions,  evident  at  50  ns  and  70  ns,  could  be  the  result  of 
nonuniformities  in  the  implosion. 
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VI.  L-SHELL  MODEL  DEVELOPMENT 


The  development  of  X-ray  dignostics  for  ions  in  the  L-shell,  using  line  emissions  from  such 
ions,  is  of  fundamental  importance  for  the  undemtanding  of  plasma  conditions.  Using  the  line 
ratio  intensities  for  two  emission  lines  arising  from  different  excited  states  of  a  given  ion  one  can 
determine  the  temperature  in  LTE.  Temperature  diagnostics  can  also  be  obtained  by  comparisions 
of  intensities  of  lines  from  neighboring  ionization  stages  of  an  element.  The  purpose  of  this  work  is 
to  develop  a  detailed  description  of  the  n=3  multiplet  levels  of  F-like  ionization  stage  of  selenium. 
We  have  already  developed  a  very  accurate  and  detailed  data  base  for  the  n=3  multiplet  levels  of 
Ne-like  selenium.  Similar  to  using  the  H  to  He  line  ratios  to  obtain  the  Ka  diagnostics,  the  F  to  Ne 
line  ratios  can  be  used  to  determine  the  L-shell  X-ray  diagnostics.  Previous  theoretical  predictions 
and  experimental  observations  indicate  substantial  F-like  population  for  selenium  plasma  and  the 
experimental  observation  of  strong  3-2  L-shell  spectra  for  this  ion  provides  the  motivation  for  a 
thorough  theoretical  investigation  of  the  F-like  selenium  ionization  stage.  A  detailed  description  of 
the  multiplet  levels  for  an  ion  such  as  F-like  selenium  requires  an  enormous  amount  of  atomic  data. 
In  order  to  accurately  determine  the  populations  of  the  atomic  states  of  an  ion  in  the  steady  state 
for  a  given  plasma  density  and  temperature,  and  to  make  comparisions  between  observed  spectra 
and  theoretical  calculations  one  needs  to  obtain  accurately  calculated  atomic  rates  populating  and 
depopulating  these  states.  Theoretical  modeling  of  the  extreme  UV  specturm  for  F-like  selenium 
in  the  experiments  at  the  Lawemce  Livermore  National  Laboratory  is  at  discrepany  with  the 
experimental  observations.  Also,  theoretical  calculations  using  plasma  hydrodynamics  simulations 
predict  strong  3-3  lines  for  F-like  ions,  however,  no  such  lines  were  observed  in  the  experiments. 
These  facts  necessitate  further  detailed  and  accurate  modeling  of  the  kinetics  of  F-like  ioniztion 
stage.  The  rates  of  the  various  atomic  processes  populating  the  states  depend  on  the  rate  coefficients 
which  are  obtained  by  calculating  the  atomic  cross  sections  and  averaging  them  over  the  electron 
distributions  for  the  appropriate  plasma  temperatures.  The  cross  sections  that  are  needed  for  the 
atomic  processes  include  collisional  excitation,  deexcitation,  ionization,  photoionization,  as  well 
as  radiative  and  dielectronic  recombination.  Under  coronal  equilibrium  conditions,  one  needs  to 
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take  into  account  recombination  processes  as  well  as  radiative  cascade  from  upper  levels  to  obtain 
various  level  populations  by  solving  a  multitude  of  different  rate  equaions.  The  level  of  the  detailed 
calculations  of  the  cross  sections  of  these  processes  as  well  as  the  determination  of  the  energy 
levels  of  the  atomic  states  involved  depend  on  the  atomic  model  required  to  accurately  represent 
the  kinetic  properties  of  the  plasma.  For  example,  in  the  ionization  balance  calculations  it  may  be 
sufficient  to  calculate  configuration  averaged  atomic  data  whereas  detailed  calculations  at  the  fine 
structure  levels  are  required  in  predicting  spectra. 

We  have  developed  the  n=3  multiplet  levels  of  F-like  ionization  stage  of  selnium  in  a  similar 
fashion  as  was  done  for  Ne~like  selenium.  Our  compct  model  contains  the  three  n=2  low  lying 
states  and  fifty  seven  multiplet  states  including  8  2p*3s,  21  2p*3p,  28  2p*3d,  three  iimer-shell 
excited  configurations  of  2p*3a,  2p*3p,  2p^3d  and  the  two  lumped  n=4  excited  state  levels.  Figure 
1  shows  this  basic  model. 

The  atomic  structure  and  cross  sections  for  this  work  were  calculated  using  various  computer 
programs  such  as  the  computer  codes  developed  at  Los  Alamos  National  Laboratory,  the  atomic  code 
of  R.  D.  Cowan^  and  the  relativistic  distorted  wave  codes  developed  by  Sampson  and  2iang.’  The 
eneigy  levels,  oscillator  strengths  and  plane  wave  Bom  (PWB)  collision  strengths  were  calculated 
using  CATS,’  a  version  of  the  atomic  structure  code  of  Cowan,  modified  by  the  group  at  Los  Alamos. 
Even  though  our  basic  atomic  model  contains  the  inner-  shell  excited  states  at  the  configuration 
level,  the  multiconfigurational  structure  calculation  includes  the  fine  structure  levels  of  these  states 
with  a  total  of  1 13  states  for  all  L-  and  M-shell  ground  and  excited  states  of  F-like  selenium.  Our 
energy  levels  compare  very  well  with  the  fully  relativisic  energies  calculated  by  Hagelstein.^  We 
have  neglected  any  configuration  (Cl)  among  higlnr  levels  and  doubly  excited  levels.  In  order  to 
check  the  effect  of  O  on  the  energy  levels,  Hagelstein  included  1622  states  which  included  313V 
doubly  excited  levels  and  found  on  significant  effect  In  Table  I  we  show  the  designations  and  the 
eneigies  of  the  levels  that  are  included  in  our  model. 

We  have  also  calculated  the  wavelengths  and  absorption  oscillator  strengths  for  some  of  the 
strong  ns2-3  transitions.  Except  for  a  few  transitions  our  absorption  oscillator  strengths  agree 
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very  well  with  those  of  Hagelstein.  When  compared  with  the  experimental  data  our  oscillator 
strengths  are  also  in  good  agreement  for  many  of  the  transitions.  Table  n  shows  a  comparision  of 
our  oscillator  strengths  with  other  data. 

From  our  previous  experience,  during  the  development  of  a  detailed  atomic  data  base  for  Ne- 
like  selenium,  we  realized  that  the  existing  Branch  codes  for  calculating  the  necessary  collision 
strengths  and  photoionization  cross  sections  for  the  model  were  quite  inadequate.  This  made  it 
absolutely  necessary  to  obtain  and  use  other  available  computer  codes  such  as  the  CATS  code  and 
the  distorted  wave  collision  codes  developed  by  Sampson  and  Zhang  and  capable  of  producing  much 
,  more  accurate  and  consistent  atomic  data  for  collision  and  photoionization.  In  order  to  compare  the 
data  obtained  using  the  in-house  capabilities  with  these  more  sophisticated  and  detailed  codes,  we 
calculated  a  few  collision  strengths  for  F-like  ions.  Since  the  branch  code  was  developed  for  K-shell 
modeling  and  the  data  for  the  L-shell  were  obtained  inconsistently,  there  were  accidental  agreemnt 
between  our  collision  strengths  and  those  obtained  by  the  other  codes  for  certain  electron  energies. 
We  thus  could  not  rely  on  the  existing  branch  codes  and  needed  to  use  these  more  developed  and 
reliable  codes  for  determining  an  accurate  level  populations  for  the  n=3  levels  of  F-like  selenium. 
The  collision  strengths  for  excitation  from  the  low  lying  n=2  to  all  the  excited  n=3  fine  structure 
levels  were  calculated  in  the  PWB  method  using  CATS.  The  results  obtained  agreed  with  other 
published  collision  stengths  except  at  very  low  energies.  The  PWB  cross  sections  for  excitation  are 
not  sufficiently  accurate  for  spectroscopic  purposes  and  specifically  near  threshold.  We  therefore 
used  the  ACE’  (Anothed  Collisional  Excittion)  code  (also  obtained  from  Los  Alamos)  which  is 
a  modified  version  of  Mann’s*  distorted  wave  (DW)  method  and  gives  more  accurate  collision 
strengths  than  the  PWB  method.  However,  since  the  DW  calculations  are  much  slower  than  the 
PWB  calculations  and  we  need  to  generate  a  large  number  of  collision  strengths,  only  the  very 
important  transitions  were  calculated  using  ACE.  ACE  uses  the  atomic  structure  data  generated  by 
CATS.  We  have  obtained  the  collision  strengths  for  excitation  from  the  n=:2  to  n=3  states  as  well  as 
for  rollisinnal  coupling  between  the  ns3  excited  states.  We  have  obtained  the  collision  strengths 

for  excitation  at  the  fine  structure  levels.  However,  these  quantities  can  be  obtained  directly  from 
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CATS  at  the  configuration  averaged  or  at  the  multiplet  levels  as  seem  necessary.  Table  III  shows  the 
collision  strengths  for  a  few  of  the  strong  2-3  transitions  in  F-like  Se  obtained  using  the  ACE  code 
in  the  DW  q)proximation.  We  also  compare  our  results  with  those  published  data  of  Sampson  and 
Zhang.  Our  DW  collision  strengths  for  the  strong  coUisional  excitations  are  within  10%  of  those 
relativistic  DW  collision  strengths  obtained  by  Sampson  and  Zhang.  The  DW  collision  excitation 
cross  sections  <rif(E)  were  obtained  from  the  collision  strengths  O,/  by  using  the  formula 


ffy  = 


ITOq*  1 

2Ji  +  lp» 


CktiE) 


where  i  and  /  refer  to  initial  and  final  states,  oo  is  the  Bohr  radius,  2Ji  •<- 1  is  the  statistical  weight  of 
the  initial  state,  £  is  the  incident  electron  energy  and  p  is  the  momentum  of  the  incident  electron. 
Our  DW  collision  excitation  cross  sections  compare  quite  well  with  those  of  Hagelstein  in  some 
cases,  whereas  in  some  other  cases  the  agreement  is  quite  poor.  This  is  shown  in  Table  IV.  It  was 
found  that  the  results  obtained  by  Hagelstein  for  the  transitions  (1-31)  and  (2-38)  which  were  at 
great  discrepancy  with  our  cross  sections  were  erroneous. 

As  already  mentioned,  we  have  also  calculated  the  coUisional  excitations  between  the  n=3 
excited  levels  and  found  that  the  2p*3d  levels  are  strongly  coupled  to  the  2p*3p  levels.  Thus 
inclusions  of  these  mcitation  rates  wiU  have  major  influence  in  the  calculations  of  population  of 
the  n=3  multiplet  levels. 

The  dielectronic  recombination  rates  from  the  low-lying  0-like  ground  levels  to  aU  the  singly 
excited  n=3  F-like  states  were  calculated  in  detail  by  using  the  HFR  (Hartree-Fock  with  relativistic 
corrections)  method  of  Cowan.  In  order  to  calculate  the  rates  for  dielectronic  recombination  (DR) 
process,  where  a  free  electron  is  captured  by  an  ion  forming  a  doubly  excited  state  which  stabilizes 
either  by  radiative  decay  to  the  singly  excited  state  of  the  next  ion  or  by  autoionizaion,  one  needs 
to  obtain  a  huge  number  of  autoionization  and  radiative  data.  For  each  doubly  excited  F-like 
levels,  we  thus  calculated  atomic  data  which  include  rates  for  aU  possible  autoionization  channels 
inriiiHing  autoionization  to  excited  states  of  the  initial  ion  and  radiative  data  such  as  osciUator 
strengths  and  radiative  decay  rates  for  aU  possible  dipole  stabilizing  transitions  from  these  doubly 
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as  well  as  singly  excited  F-like  states.  Our  DR  data  also  necessarily  include  the  energies  for  these 
large  number  of  doubly  and  singly  excited  states. 

The  DR  branching  ratios  from  the  low  lying  n=2  0-like  levels  to  each  specific  fine  structure 
levels  of  the  singly  excited  n=3  F-like  states  were  calculated  explicitly  for  low  lying  doubly  excited 
states.  For  the  high  lydberg  states  with  n>10,  the  autoionization  rates  were  iq)proximated  by  using 
the  1/n’  fall  off  of  these  rates  while  the  radiative  rates  remained  constants.  Our  DR  branching 
ratios  were  calculated  for  a  large  number  of  intermediate  resonance  states  in  the  single  configuration 
approximation.  Approximate  relativistic  corrections  are  made  to  both  the  radial  wave  functions  as 
well  as  to  the  singly  and  doubly  excited  states  energies.  Figure  2  shows  the  DR  rate  coefficients 
to  some  of  the  n=3  levels  of  F-like  Se.  We  have  previously  calculated  the  DR  branching  ratios  for 
several  F-like  ions  including  Ar*,  Fe^*,  Se"  and  Mo"  ions  and  obtained  scaling  relations  using 

the  DR  data  for  recombination  from  At*,  H",  Fe“  and  Se"  ions  using  three  or  four  coefficients 
polynomial  fits  to  our  data.  These  results  ate  already  in  press  and  will  be  published  in  the  ‘Atomic 
Data  and  Nuclear  Data  Thbles*.^  The  DR  data  for  recombination  from  0-like  Mo  was  calculated  to 
investigate  the  validi^  of  extrapolating  the  Z-scaling  beyond  Se.  Except  for  a  few  data  points,  the 
Z-scaled  data  using  the  data  for  At*,  H",  Fe",  and  Se"  ion  predict  the  Mo"  data  which  are  within 
a  few  percent  of  the  explicitly  calculated  data.  This  and  the  fact  that  our  DR  rate  coefficients  for 
0-like  and  F-like  ions  agree  quite  well  with  other  delailed  and  more  complex  calculations,  give  us 
confidence  in  our  use  of  Cowan’s  HFR  method  of  calculations  of  DR  rates. 

Presently  our  L-shell  data  base  for  F-like  selenium  does  not  include  any  radiative  recombination 
rates.  We  will  in  the  near  future  upgrarte  our  calculations  by  including  such  rates.  We  also  intend  to 
upgrade  our  ionization  rates  for  ionization  from  Ne-like  to  F-like  selenium  very  soon.  We  already 
have  the  computer  codes  (one  from  Perm  state  by  Sampson  and  Zhang,  and  the  other  from  Los 
Alamos)  available  to  compute  improved  ionization  rates.  In  the  case  of  Ne-like  Se,  it  was  seen  that 
the  n=3  levels  of  this  ion  were  substantially  populated  by  indirect  or  resonance  excitation  (RX). 
We  therefore  plan  to  compute  such  RX  rates  and  include  them  to  the  direct  coUisional  excitation 
rates  for  excitation  from  ground  levels  to  dw  nss3  levels  of  F-like  selenium. 
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Table  I.  State  configurations,  total  angular  momenta  J  and  energies  relative  to  the  ground  state  for  the  lowest  65  states 
of  F-like  selenium.  Notation  a[b]  means  ax  10^. 


level 

jj-state 

J 

Energy  (eV) 

level 

jj-state 

J 

Energy  (eV) 

1 

[2s»2p*„2p5„ls/a 

3/2 

0 

33 

t2s>2p>,,2p»„la3d3/a 

1/2 

1.616[3] 

2 

[2s»2p*/,2pS/,]i/a 

1/2 

4328[1] 

34 

[2s»2p»,22p»/,]a3d5/a 

9/2 

1.618[3] 

3 

I2s2p»/,2pJ„li/a 

1/2 

2.131[2] 

35 

[2s»2p;„2pS/,]a3d3/a 

7/2 

1.619[3] 

4 

[2s*2pJ/j2p5/,]a3s 

5/2 

1.50U3] 

36 

[2s>2pJ,22pS/2]a3d5/a 

1/2 

1.623[31 

5 

(2s^2p»/,2p5/,]a3s 

3/2 

1.506[3] 

37 

[2s»2pJ/22p»/3]a3d5/a 

3/2 

1.627[3] 

6 

I2s»2pJ„2p»/,]o3s 

1/2 

1.519t3] 

38 

[2s*2pi/a2p5^j]a3p3/3 

1/2 

1.628[3] 

7 

[2s*2pj/a2l^yalt 

3A1 

1J42[3] 

39 

[2s»2p»/22p|/2]33d5/3 

5/2 

1.63013] 

8 

[2s*2pi/a2i^/,li3s 

1/2 

1.546P1 

40 

[2s»2p»,22i^/2lo3d3/a 

3/2 

1.635[3] 

9 

I2s»2pj„2p»,,]a3pi/a 

3/2 

1.547[3] 

41 

[2s»2p?,22p5,2lo3d5/a 

5/2 

1.638[3] 

10 

[2s^2pJyj2pj^2]a3pi/a 

5/2 

1.548[3] 

42 

[2s^2pi/32p3^2]i3d3/a 

1/2 

1.65313] 

11 

[2s’2pi/a2i^/,la3s 

5/2 

1.555[3] 

43 

I2s*2p5^2lo3pi/3 

1/2 

1.656[3] 

12 

[2s»2pi/a2p>„]a3s 

3Al 

1.556[3] 

44 

[2s’2pi/a2p’^2]i3d3/a 

3/2 

1.657[3] 

13 

[2s*2pJ/j2i^/j]a3p3/a 

1/2 

1357[3] 

45 

t2s*2pi/a2p5y2lx3d5/a 

7/2 

1.658[3] 

14 

[2s»2pJ/22p»„]a3p3/a 

5/2 

I.558i3] 

46 

[2s’2pJ/2Jo3ps/a 

3/2 

1.659[3] 

15 

[2s*2pJ,j2i^/2la3pa/a 

7/2 

1,559[3] 

47 

[2s*2pi/a2p|/2]i3d5/a 

5/2 

1.661[3] 

16 

I2s*2Pi/a2|^/a3o3pi/a 

1/2 

1J66[3] 

48 

[2s’2pi/a2||^/2li^9/3 

5/2 

1.662[3] 

17 

[2s*2p5/22p5/,]a3ps/a 

3/2 

lJ72t3] 

49 

[2s’2pi/a2p^^2]i^‘^>/3 

3/2 

1.66313] 

18 

f2s*2p»/,2r^„]o3p3/a 

3/2 

1377t3] 

50 

[2s’2pi/a2i^^,]a3d3/a 

7/2 

1.668[3] 

19 

[2s*2pi/a2p|,,]i3pi/a 

1/2 

1J86[3] 

51 

[2s»2pa/a2i^/2]»3d»/* 

9/2 

1.670[3] 

20 

[2s’2pi/a2l^^alt3pi/a 

3/2 

1.591  [3] 

52 

[2s’2pi/a2p3^2]3^B/3 

5/2 

1.672[3] 

21 

t2s*2pi/a2i^,3]i3p3/a 

5/2 

1.598[3] 

53 

[2s’2pi/32p|^3]a3d3/3 

1/2 

1.672[3] 

22 

I2s*2pi/a2p5/j]a3pi/a 

5/2 

1.600[3] 

54 

[2s*2pi/32p5/2l*3d5/3 

7/2 

1.675[3] 

23 

[2s’2pi/a2i^^,]i3ps/a 

1/2 

1.601  [3] 

55 

[2s*2pj/32p|/2]23d3/3 

3/2 

1.67S[3] 

24 

[2s’2pi/a2p3/alx3p3/a 

3/2 

1.602[3] 

56 

I2s*2pi/32p3/2lB3d3/3 

5/2 

1.677[3] 

25 

I2s>2p</,]o3s 

1/2 

1.603[3] 

57 

[2s*2pi/32i^/2l»3d5/a 

3/2 

1.684[3] 

26 

[2s’2pi/a2i^^,]a3p3/a 

7/2 

1.61 1[3] 

58 

[2s^2pa/32p3^3]a3d5/a 

1/2 

1.686[3] 

27 

t2s*2pt/a2i^^jla3p3/a 

3/2 

1.611[3] 

59 

I2s»2pJ,2]o3d8/3 

5/2 

1.719[3] 

28 

[2s»2p»/,2[4„la3d,/a 

5/2 

1.614[3] 

60 

[2s»2pS/2]o3d3/3 

3/2 

1.72313] 

29 

[2s»2p5/22p|/2]a3d3/a 

3/2 

1.614[3] 

61 

2s2p*3s 

1.705[3] 

30 

I2s*2pi/a2p|/j]33p3/a 

5/2 

1.615[3I 

62 

2s2p*3p 

I.758f3] 

31 

[2s’2pi/a2p3^,]a3pi/a 

3/2 

1.615[3] 

63 

2s2p®3d 

1.820[3] 

32 

[2s*2pJ/j2pJ/3]a3d5/a 

7/2 

1.615[3] 

64 

65 

2s*?pMl 

2s2p®41 

2.102[3] 

2.282[3] 
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Table  n.  Wavelengths  and  absorption  oscillator  strengths  for  strong  2p-3s  and  2p-3d  transitions  in  F-like  Se. 


Transitions 

Ji-Jr 

A  (A) 

(this  work) 

HRS 

f 

(this  work) 

f 

(Ref.  3) 

2-8 

1/2-1/2 

8.252 

8.253 

1-5 

3/2-3/2 

8.234 

8.237 

2-12 

1/2-3/2 

8.195 

8.196 

0.094 

0.089 

1-11 

3/2-5/2 

7.974 

7.978 

0.054 

0.049 

2-25 

1/2-1/2 

7.950 

7.948 

0.044 

0.089 

1-37 

3/2-3/2 

7.621 

7.625 

0.181 

0.185 

1-39 

3/2-5/2 

7.606 

7.610 

0.382 

0.381 

1-41 

312-Sn 

7J69 

7.570 

0.268 

0.243 

2-57 

m-zn 

7.558 

7.558 

0.977 

0.984 

2-58 

1/2-1/2 

7.547 

7346 

0.774 

0.770 

1-48 

za-sn. 

7.458 

7.464 

0.181 

0.205 

1-53 

zn-m 

7.414 

7.419 

0.227 

0.220 

1-55 

zn-zn 

7.401 

7.405 

0.4% 

0.497 

1-56 

3/2-5/2 

7.395 

7.398 

0.792 

0.773 

2-64 

1/2-3/2 

7.360 

7.380 

1.019 

0.979 

1-57 

zn-zn 

7.364 

7.366 

0.145 

0.130 
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lUile  m.  Collision  strengths  for  excitation  from  some  of  the  strong  2-3  transitions  in  F-like  selenium.  Notation  a[b] 
means  axlO^. 


Transitior 

1  Transition  Energy(eV)  E/Eth 

lUsWork 

Ref.  7 

1-17 

1.572(3] 

1.058 

1.60(-2] 

1.46(-2] 

1.260 

1.63(-2] 

1.50(-2] 

1.646 

1.68(-2] 

1.55(-2] 

2.360 

1.75(-2] 

1.62(-2] 

3.650 

1.81(-2] 

1.70(-2] 

5.849 

1.87(-2] 

1.78(-2] 

1-31 

1.615(3] 

1.050 

5.47(-2] 

5.77(-2] 

1.250 

5.62(-2] 

5.95(-2] 

1.629 

5.84(-2] 

630(-2] 

2320 

6.06(-2] 

6.48(-2] 

3.58 

635(-2] 

6.75(-2] 

5.722 

6.37(-2] 

6.98(-2] 

2-38 

1J85(3] 

1.051 

2.17(-2] 

1.92(-2] 

1360 

233(-2] 

1.97(-2] 

1.640 

2.31  (-2] 

2.05(-2] 

2350 

239(-2] 

2.14(-2] 

3.630 

2.46(-2] 

233(-2] 

5.802 

230(-2] 

230(-2] 

2-43 

1.613(3] 

1.050 

2.!6(-2] 

2.19(-2] 

1350 

230(-2] 

235(-2] 

1.630 

237(-2] 

232(-2] 

2320 

234(-2] 

2.41(-2] 

3.630 

2.40{-2] 

230(-2] 

5.726 

2.43(-2] 

238(-2] 

1-56 

1.677(3] 

1.050 

5.99[-2] 

6.12(-2] 

1340 

6.81(-2] 

6.97(-2] 

1.610 

837(-2] 

8.44(-2] 

2370 

1.05(-1] 

1.07(-1] 

3.490 

1.37(-1] 

1.40(-1] 

5350 

1.74(-1] 

1.80(-1] 

2-57 

1.640(3] 

1.050 

3.88(-2] 

3.97(-2] 

1350 

4.34(-2] 

434(-2] 

1.620 

537(-2] 

5.49(-2] 

2.300 

6.81(-2] 

6.98(-2] 

3340 

8.83(-2] 

9.10(-2] 

5.650 

1.12(-1] 

1.17(-1] 

2-64 

1.680(3] 

1.050 

3.98(-2] 

3.95(-2] 

1.240 

4.17(-2] 

4.46(-2] 

1.600 

5.34(-2] 

534(-2] 

2.270 

6.78(-2] 

6.74(-2] 

3.480 

8.77(-2] 

8.77(-2] 

5340 

l.lU-1] 

1.13[-2] 
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Table  FV.  Collision  excitation  cross  sections  (cm’)  for  excitation  from  some  of  the  strong  2-3  transitions  in  F-like 
selenium.  Notation  a{b]  means  axIC*^ 


Transition 

Transition  Energy(eV) 

E/Etfc 

TTiisWork 

Ref  .  3 

1-17 

1.572(3] 

1.13 

2.711(-21] 

13191-21] 

1.32 

2.365(-21] 

1.3401-21] 

1.64 

1.950(-21] 

1.1251-21] 

221 

1.460(-21] 

8.5971-22] 

3M 

9.730(-22] 

5.9171-22] 

1-31 

1.615(3] 

1.13 

9.080(-21] 

1.7311-20] 

IJl 

8.020(-21] 

13231-20] 

1.63 

6.630(-21] 

13671-20] 

2.32 

4.840(-21] 

9.4801-21] 

3.48 

3.320(-21] 

63731-21] 

2-38 

1.585(3] 

1.13 

7320(-21] 

1.9041-21] 

132 

6.440(-21] 

1.7121-21] 

1.64 

5320(-21] 

1.4661-21] 

236 

3.980(-21] 

1.1431-21] 

332 

2.640(-21] 

7.9661-22] 

2^3 

1.613(3] 

1.12 

7.190(-21] 

43951-21] 

131 

6390(-21] 

3.8031-21] 

1.62 

5300(-21] 

3.1971-21] 

234 

3.860(-21] 

2.4261-21] 

3.48 

23501-211 

1.6411-21] 

1-56 

1.677(3] 

1.12 

1.000{-20] 

1.0401-20] 

130 

9.680(-21] 

1.0111-20] 

1.61 

9.160(-21] 

9.6281-20] 

2.19 

8.310[-21] 

8.7821-21] 

338 

7.0801-21] 

73381-21] 

2-57 

1.640(3] 

1.12 

13301-20] 

1.3601-20] 

130 

1.2801-20] 

1.3171-20] 

1.62 

13101-20] 

13491-20] 

232 

1.0901-20] 

1.1291-20] 

3.44 

93201-21] 

93821-21] 

2-64 

1.680(3] 

1.12 

1.3301-20] 

1.3021-20] 

130 

13701-20] 

13531-20] 

1.60 

1.1901-20] 

1.1831-20] 

2.19 

1.0801-20] 

1.0681-20] 

338 

9.1081-21] 

9.0901-21] 
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